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INTRODUCTION 


The  primary  objective  of  this  effort  has  been  the  investigation  of  dominant  and 
higher-order  model  characteristics  of  microstrip  lines  and  fin  lines  that  find  application  in 
millimeter-wave  integrated  circuits.  The  investigation  of  higher-order  modes  is  important 
for  the  analysis  of  discontinuity  problems  in  planar  transmission  lines:  these  problems  have 
been  the  focus  of  attention  in  this  project  during  the  entire  grant  period.  The  study  of  the 
discontinuity  problems  is  of  interest  from  the  point  of  view  of  designing  practical 
components,  e.g..  filters,  impedance  matching  networks,  and  transitions  between  rectangular 
waveguides  and  planar  transmission  lines.  The  higher-order  modes  have  been  analyzed 
using  the  spectral  domain  technique  introduced  by  Iloh  and  Mittra  [l.  2].  The  propagation 
constant,  characteristic  impedance  and  field  configuration  are  obtained  for  the  dominant  and 
higher-order  modes  in  planar  lines.  The  model  characteristics  are  then  employed  to  analyze 
some  typical  discontinuities  in  planar  lines  using  the  mode-matching  procedure  [3.  4], 
which  has  been  extensively  employed  in  the  past  for  the  solution  of  similar  problems  in 
other  types  of  waveguides.  However,  the  situation  with  regard  to  the  application  of  the 
mode-matching  procedure  is  found  to  be  quite  different  for  the  case  of  planar  waveguides 
when  compared  to  that  of  rectangular  guides  because  accurate  generation  of  higher-order 
model  fields  becomes  a  formidable  problem  in  planar  guides,  whereas  this  task  is  quite 
straightforward  for  rectangular  guides.  Thus,  it  becomes  important  to  investigate 
effectiveness  of  the  mode-matching  procedure  when  only  a  relatively  small  number  of 
modes  are  available  and  the  higher-order  modes  are  known  only  approximately. 

In  an  attempt  to  enhance  the  accuracy  of  the  higher-order  mode  computation,  the 
singular  integral  equation  technique  is  used.  This  technique  was  originally  employed  by 
Mittra  and  Iloh  [5]  for  the  analysis  of  planar  transmission  lines.  This  technique  is  further 
developed  in  this  work  and  the  results  obtained  via  the  application  of  the  singular  integral 
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equation  technique  are  compared  with  those  derived  from  the  spectral  domain  approach. 
Another  interesting  problem  in  planar  transmission  lines,  viz.,  the  coupled  multiconductor 
lines  is  analyzed  in  this  work  using  the  spectral  domain  approach. 

To  the  best  of  our  knowledge  a  full-wave  coupled  line  analysis  of  n  planar 
transmission  lines  has  not  appeared  elsewhere  in  the  literature. 


CHAPTER  2 


UNIFORM  MICROSTRIP  ANALYSIS 


2.1  Introduction 

The  analysis  of  various  printed  circuits  has  been  of  interest  for  a  number  of  years. 
We  find  it  of  interest  here,  because  the  first  step  of  the  mode  matching  procedure  is  to 
generate  the  dominant  and  higher-order  modes  in  a  uniform  microstrip.  A  cross  section  of  a 
shielded  microstrip  is  shown  in  Figure  2.1.  This  and  related  structures  have  been  analyzed 
by  a  number  of  workers  using  a  variety  of  techniques.  These  techniques  include  various 
quasi-static  methods  [6-9],  nonuniform  discretization  of  integral  equations  [10],  a  modified 
Weiner-Hopf  technique  [ll],  a  singular  integral  equation  technique  [5].  and  the  spectral 
Galerkin  technique  [1.2.12-17],  The  method  that  is  simplest  to  implement  while  giving 
excellent  results  is  the  spectral  Galerkin  technique.  This  is  the  method  considered  in  this 
chapter. 

In  order  to  allow  mode  matching  work  to  a  high  degree  of  accuracy,  it  is  necessary 
to  calculate  as  many  evanescent  modes  as  possible.  The  feasibility  of  generating  a  large 
number  of  higher-order  modes  imposes  the  largest  constraint  upon  the  accuracy  of  the 
mode  matching  solutions.  To  date,  there  is  little  information  available  on  evanescent  modes 
in  a  shielded  microstrip,  although  the  propagating  modes  have  been  well  studied. 

In  this  chapter,  we  present  data  on  the  dominant  and  first  few  evanescent  modes  of 
a  microstrip,  showing  dispersion  curves,  characteristic  impedance  calculations,  and  field 
distrbutions.  Let  us  proceed  now  with  the  analysis  of  a  uniform  microstrip. 

2.2  The  Spectral  Domain  Immitance  Approach 

In  order  to  find  information  about  waveguide  modes,  the  quantity  of  most 
immediate  interest  is  the  propagation  constant.  l3.  To  find  (3.  a  matrix  equation  must  be 


found  of  the  form 


Z21w  Z„(0)  Jz 
z'  O)  z„(fl)  /, 


£z(0) 

i,  (jb) 


where  /  and  £  are  the  Fourier  transforms  of  the  current  and  electric  field  in  the  plane  2«0. 
and  the  Z 's  form  the  components  of  a  dyadic  Green's  function.  In  order  to  find  this  Green's 
function,  we  use  the  spectral  domain  immitance  approach.  This  approach  is  a  method  of 
generating  the  dyadic  Green’s  function  in  a  straightforward  manner,  that  will  be  useful  for 
many  different  kinds  of  printed  circuits,  including  microstrip,  fin  line,  and  coplanar 
waveguides.  The  work  in  this  section  follows  Itoh  [13]. 

We  begin  by  setting  up  a  scalar  Helmholz  equation  in  each  of  the  regions  1  and  2.  as 
shown  in  Figure  1.  Assuming  anej(“  dependence,  we  obtain 


(*  ■ y ) 

(  V2  +  €.*„2  )  =  0 

<M*.y) 


(2.2a) 


where 


eo  Region  1 

€‘  —  e,  Region  2 


(2.2b) 


As  the  equations  now  stand,  the  TE}  and  TMy  fields  are  coupled.  In  order  to  decouple  these 


equations  we  must  work  in  the  Fourier  transform  domain.  Thus,  for  example. 


4>  (n.  .y  )  =  f  <t>  U  .y  )  e 1  a" '  dx 


0  (*  -y  )  "  2b  ^  e~ian* 


_  (n  -  Vi  )  IT 

*"  b 


(2.3a) 


(2.3b) 


(2.3c) 


The  above  choice  of  a,  is  suitable  for  modes  that  are  even  in  Jz .  such  as  the  dominant 
mode.  This  choice  enforces  the  boundary  conditions  at  x  =  ±b  .  For  modes  that  are  odd  in 
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Jz  we  choose  a,  =  ntr/b  .  In  the  transform  domain,  the  Helmholtz  equations  become 


,2  ) 

(-a„2+iT-)JJ  +  «,it2)  -  =° 

3y  *(«) 


o  6  (n  ) 

(A-y?)  -  =0 

3>’2  *(n) 


where 

V,2  =  a2  +  02-e,*o2  (2.6) 

Thus,  in  the  transform  domain  we  can  reduce  the  problem  to  two  transmission  line 
equations. 

Let  us  now  introduce  the  transformations 


2 

jV.  N, 

V 

a: 

=  Nx  -Nt 

u 

where 


Nz  =  — - - -  =  sin  0 

VaJ+02 


(2.8a) 


N,  ~  —  a  =  COS  0 

V  <*n  +  P2 


(2.8b) 


This  transformation  is  shown  is  Figure  2.2.  With  this  transformation,  we  see  that  all  scalar 
potentials  and  fields  are  of  the  form 

QO 

l  (n  ,v  )  e  ■  f <t>(x  ,v  )  e  '  +  ’  (2.9) 


The  transformed  equations  are  dependent  only  upon  v;  hence  -JL  =  0.  This  decouples 

0“ 


! 

I 

i 

> 

[  8 

i 

> 

[  Equation  2.5  into  TEy  components  UJy  ,  Eu  .  and  H v  )  and  TMy  components  (Ey  .  Ev .  and 

|  H „  ).  We  now  show  how  this  separation  takes  place,  and  that  the  TEy  components  are  due 

to  J u .  and  the  TMy  components  are  due  to  J  v . 

Let  us  derive  the  expressions  for  the  field  components.  To  begin,  we  express  the 
vector  potentials  as 


F ,  -y 4>i(x.y) 

(2.10a) 

Ai  =y+,ix.y) 

(2.10b) 

Adapting  some  results  from  Harrington  [18]  on  the  separation  of  TE  and  TM  fields,  and 

using  the  fact  that  —  y,  and  —  0.  we  obtain  the  TE  components 

0 

(2.11a) 

ll 

(2.11b) 

~  _  y,  di, 

jd>Ho  0V 

(2.11c) 

H-  ■  ,  {y? + k'2)*' 

(2.1  Id) 

ii 

© 

(2.1 te) 

and  the  TM  components 

r  _ 

v  /<**,  dv 

(2.12a) 

£>7^r(v'+M)i 

(2.12b) 

Eu  =0 

(2.12c) 

From  these  equations,  we  can  set  up  equivalent  transmission  lines  for  the  TE  and  TM  cases. 
The  characteristic  impedance  of  each  line  is 


11 

where 

Z~«  =  N*Z'  +  N*Zh  (2.17a) 

ilx  =  Z‘«  =  Nx  Nt  (  -  1’  +  ih  )  (2.17b) 

Z'„  =  NX2Z'  +  Nt2Zh  (2.17c) 

This  is  the  equation  that  must  be  solved  to  find  0.  Once  0  is  found,  all  other  quantities  of 
interest,  such  as  field  configurations  and  characteristic  impedance,  may  be  found. 

It  is  helpful  at  this  point  to  briefly  discuss  the  applicability  of  the  spectral  domain 
immitance  approach  to  structures  related  to  the  microstrip.  Since  the  equations  are 
separable  in  the  transform  domain,  we  see  that  if  we  have  multiple  layers  of  dielectric,  this 
will  be  equivalent  to  having  extra  sections  in  the  transmission  line  model.  Furthermore, 
multiple  strips  can  be  represented  as  extra  current  sources,  and  slots  can  be  represented  by 
voltage  sources.  Therefore,  we  can  use  this  method  of  generating  the  dyadic  Green's 
[  function  for  a  large  class  of  planar  waveguide  structures.  These  include  fin  line,  coplanar 

waveguide,  and  slot  line.  The  generality  in  this  procedure  for  generating  the  Green  function 
maintains  the  generality  of  the  mode  matching  procedure. 

I 

23  The  Spectral  Galerkin  Technique 

Now  that  we  have  found  Equation  2.16,  we  must  solve  it  for  0  with  the  spectral 
Galerkin  technique.  The  work  in  this  section  follows  Schmidt  and  Itoh  [15]. 

We  proceed  by  expanding  the  currents  on  the  strip  as  a  sum  of  basis  functions 

;  s 

■M*)*  Z  c,£,(x)  (2.18a) 

1=1 

y 

Jx  (*)“  (2.18b) 

’  1*1 


These  basis  functions  are  now  Fourier  transformed  according  to  Equation  2.3  and 
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substituted  into  Equation  2.16.  By  taking  the  inner  product  of  the  resulting  equation  with 
(,  ( x  )  and  Tj,  (x  ).  we  can  eliminate  the  E-fields  in  this  equation.  This  is  true  because  ( x  ) 
and  tj,  (x  )  are  non-zero  only  where  Et  and  £,  are  constrained  to  be  zero,  since  they  are 
tangential  to  the  strip.  Since  this  is  true  in  the  space  domain.  Parseval's  theorem  [  1 8] 
guarantees  that  it  will  also  be  true  in  the  spectral  domain.  Hence,  we  obtain  the  equations 


where 


C,  +  =0  p=\....M 

1=1  1  =1 


=0  q=  1 . N 

i*i  j*i 


K?,  =  I  &(n)z‘„(n.0)£(n) 

n  =— oo 


KPJ  =  Z  (p(n)  Z:i(n  .0)  T)j(n) 


AT”  =  Z  Vn)z'xz(n.0)£,(n) 


=  z  (n  )  (n  .0)  (n  ) 


(2.19a) 


(2.19b) 


(2.20a) 


(2.20b) 


(2.20c) 


(2.20d) 


A  solution  to  this  equation  exists  if  and  only  if  the  determinant  of  the  coefficient  matrix  is 
equal  to  zero.  Thus. 


o  —  f  (0)  —  KCC 
/  Kxz  Kxx 


(2.21) 


where,  for  example.  K::  is  the  matrix  of  coefficients  given  in  Equation  2.20a.  We  may  now 
find  the  zeros  of  f(0)  with  a  zero-finding  algorithm  such  as  Newton's  method,  thus 
completing  the  solution  for  0. 


N  .•  *.  \  N  *»  ,N  *. 


V  N  S.  *W  N.  *•- 
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2.4  Basis  Functions 

One  of  the  most  important  aspects  of  the  spectral  Galerkin  procedure  involves  the 
choice  of  basis  functions.  A  good  set  of  basis  functions  can  greatly  increase  the  efficiency  of 
the  solution,  as  well  as  its  accuracy. 

The  basis  functions  should  satisfy  three  conditions.  First,  they  must  be 
analytically  Fourier  transformable,  since  they  are  used  in  the  transform  domain  in 
Equation  2.16.  Second,  they  should  have  a  shape  that  is  similar  to  that  for  the  current  we 
expect  to  find  on  the  strip.  Having  this  property  will  reduce  the  size  of  the  matrix 
equation.  Finally,  they  should  have  the  properties  that  Jt(x)  haS  a  1/VjT  singularity  at 
the  strip  edges  and  that  Jx  (x  )  is  zero  at  the  strip  edges.  This  is  the  result  expected  for  the 
current  parallel  and  perpendicular  to  a  knife  edge.  By  not  satisfying  this  last  condition, 
more  basis  functions  will  be  required  to  represent  the  currents,  resulting  again  in  an 
increased  matrix  size. 

One  appropriate  choice  of  basis  functions  for  modes  that  are  even  in  Jt(x)  is 
suggested  in  [  15]  as 


£  (x)~  cos  [  (t  -1)  v  (x/j  +  1)} 
Vl  —  (x  A)2 

,  \  sin  [i  it  (x  Is  +1)] 

T),  U  1  -  - ===== - 

Vl  -ix/s)2 


(2.22a) 


(2.22b) 


These  functions  satisfy  all  the  above  criteria.  In  particular,  note  that  (x  )  has  a  1/x 
singularity  at  x-±s  and  that  7},  (x  =  ±r  )  =  0.  which  can  be  shown  by  l'Hopital's  rule. 
The  first  few  functions  are  shown  in  Figures  2.4  and  2.5.  The  Fourier  transforms  of  these 
f unctions  are 

f  (n  )  =  [7,(a„  s  +  U-l)rr)  +  /0(a„r  -  (i -!)«■)]  (2.23a) 


n  (n  )  -  - — l?  S  n  [Jr,  ( a„s  +  i  n)  -  J„( a,  s  -  i  tr)] 

zy 


(2.23b) 
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where  j  =  •J—'. 1  and  Ja  (x  )  is  the  zeroth  order  Bessel  function.  A  formula  that  is  useful  in 


deriving  the  above  transforms  comes  from  the  Bateman  Manuscripts  [19].  Hence. 
“  T  (x  )  e  ~ixy 

f  *  .  dx  =  (-1Y  in  ttJ„  ( y  ) 

4,  Vi_x2 


(2.24) 


where  T„(x)  is  the  nth  order  Chebychev  polynomial  in  the  interval  lx  I  <  1  and  zero 
elsewhere,  and  Ta  (x  )  =  1. 


2-5  Characteristic  Impedance 

The  characteristic  impedance  of  a  microstrip  is  a  useful  parameter  to  calculate  for 
several  reasons.  It  serves  not  only  as  a  useful  design  tool,  but  also  as  a  check  on  the 
accuracy  of  the  dominant  mode  propagation  constant.  In  addition,  the  inner  product  we  use 
in  the  mode  matching  procedure  uses  an  inner  product  that  is  very  similar  to  characteristic 
impedance.  Hence,  it  will  be  useful  to  compare  our  calculated  values  of  characteristic 
impedance  to  previous  calculations  [11.12]. 

At  this  point,  mention  should  be  made  of  the  ambiguity  inherent  in  this 
characteristic  impedance  calculation.  The  characteristic  impedance  of  a  transmission  line  is 
usually  thought  to  be  a  characteristic  of  TEM  transmission  lines.  The  dominant  mode  of 
the  microstrip  is  only  quasi-TEM.  however,  so  there  are  several  definitions  of  characteristic 
impedance  in  this  case  that  make  sense  and  yield  similar,  but  unequal  results.  A  number  of 
authors  have  discussed  the  relative  merits  of  the  various  definitions  [20-23].  These 
definitions  include 


(2.25) 


Z 

0  / 


(2.26) 
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where 


Z. 


w 

2  P 


h-t  b 

P  =  Vi  Re  f  f  E  X  H‘  ■  z  dx  dy 

— f  “6 


(2.27) 


(2.28) 


/ 


/ 


/.  (x  ) 


(2.29) 


o 

V  =  -  f  E,(x  =0.y  )  <fy  (2.30) 

— f 

These  are  the  so-called  power-current,  voltage-current,  and  power-voltage  definitions.  It 
turns  out  that  the  definition  most  widely  accepted  is  the  power-current  definition,  and  that 
is  what  we  calculate  here.  For  comparison,  we  also  calculate  the  characteristic  impedance 
with  the  voltage-current  definition.  Thus,  we  need  to  calculate  P.  I.  and  V. 

First,  let  us  calculate  I  and  V.  Proceeding  from  Equations  2.29  and  2.30.  it  is 
straightforward  to  show  that 

f  =  Z  c.I.U  =0)  (2.31) 

I  =1 


V 


2 T-  L  jEAn.y)dy 

n  =— oo  — f 


(2.32) 


Next,  we  proceed  to  the  power  calculation.  Beginning  with  Equation  2.28.  and 
subsequently  using  Parseval's  theorem  [24],  we  obtain 

P  I  /(*)  (2.33) 

w  n  *-oo 


where 

H  —t 

/ (n  )  =  /  (£,  HI  -  E,  H[ )  dy 


(2.34) 
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The  task  now  becomes  one  of  finding  expressions  for  the  field  components  in  the  transform 


domain.  In  order  to  find  these,  we  first  need  to  solve  the  transmission  line  problem  shown 


previously  in 

Figure  2.3.  The  solution  is  straightforward,  and  can 

be  adapted  from  a 

number  of  standard  textbooks,  for  example.  Mayes  [25].  This  results 

in  the  following  set 

of  equations 

in 

.y  )  =  A,rE(n  )  Xf(y  ) 

(2.35a) 

Em  (n 

,y  )  =  ~A,rHn  )  X,s(y  ) 

(2.35b) 

Hm  ( n 

■y  )  =  A,™(n  )  X,s(y  ) 

(2.35c) 

(n 

C: 

II 

a* 

"  -i 
>» 

3" 

X 

' — ✓ 

(2.35d) 

where 

cosh  y^v  —  (h  —  t  )) 

xf  = 

sinh  yx(h  —  t )  t  =  1 

cosh  y2(y  +  t )  i  =  ■> 

sinh  y2t 

(2.36a) 

sinh  y^y  —  (h  —  i  )) 

X*  = 

sinh  y |(/i  —  r )  i  =  1 

sinh  y2(y  +  t  )  t  -  2 

sinh  y2r 

(2.36b) 

and 

■4?  a 

— Z  f£  sinh  y,(/i  —  i  ) 

(2.37a) 

A  \- 

Z  2£  sinh  y^f 

A 

= 

— Z  jA/sinh  yx(h  —  t  ) 

Z  2V/  sinh  yy 

12. J  fO) 

(2.37b) 
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Hy<  (n.y)  = 


V  On2  +  P2 

mT~ 


AFzpx* 


(2.40e) 


H „  (n  .y  )  =  |  iV2  A,rf  +  N,  A™ \  *,c  (2.40f ) 

The  v  dependence  is  contained  in  Xf  and  X /.  These  fields  can  now  be  substituted  into 
Equation  2.7  to  calculate  f(n).  In  order  to  carry  out  the  integration  of  f(n).  we  need  only 
integrate  over  the  products  of  hyperbolic  cosines  and  sines,  which  can  be  carried  out 
analytically.  This  concludes  the  power  calculations  and.  hence,  the  characteristic 
impedance  calculation. 


2.6  Fin  Line  Calculations 

Although  it  is  very  easy  to  verify  the  dominant  mode  microstrip  calculations,  very 
little  data  exist  for  the  evanescent  modes.  We  would,  however,  like  to  have  a  comparison 
for  our  calculated  evanescent  modes.  It  turns  out  that  fin  line  evanescent  modes  have  been 
calculated.  In  order  to  verify  our  evanescent  modes,  therefore,  the  best  we  can  do  is  to 
alter  our  program  to  calculate  fin  line  evanescent  modes,  and  compare  these  to  those  in  the 
literature  [17],  A  diagram  of  fin  line  is  shown  in  Figure  2.6. 

In  order  to  adapt  the  spectral  domain  immitance  approach  to  the  fin  line,  we  need  to 
"se  a  different  matrix  equation.  This  is  of  the  form 


Ex 

J, 

ylt 

£ 

i i 

(2.41) 


where  E ,  and  E.  in  the  fin  line  are  analogous  to  Jz  and  J x  in  the  microstrip,  respectively. 
In  order  to  find  the  admittances,  we  must  set  up  transmission  lines  analogous  to  those  in 
Figure  2.3.  These  are  shown  in  Figure  2.7.  From  these  transmission  lines  we  obtain  the 
dyadic  admittance  matrix 


Y 


KX 


=  N:  Y’  +  N?  Yh 


'i 


W 


(2.42a) 


Y„  =  Nt  -Y'  + 


Nz2  Y‘  +  Nx2  Y* 


2 1  lanh  y\\h  —a  —l  ) 

1  Z?1  +  2 [v;  tanh  y^  tanh  y2t 


tanh  Via  +  Z 


lanh  y\{h  —a  —t 

1  Z%E  +Zf£  tanh  yta  lanh  y^ 
ZF  Z^  tanh  yqa  +  ZiE  tanh  y2f 


(2.43b) 


The  matrix  equation  is  now  solved  in  a  manner  very  similar  to  the  microstrip  matrix 
equation.  We  expand  Ex  and  E.  in  the  same  basis  functions  used  previously  for  J t  and 
J ,  .  respectively.  These  basis  functions  were  given  in  Equation  2.22. 

Finally,  we  must  reconsider  the  definition  of  the  Fourier  transform  before  we  solve 
the  fin  line  problem.  In  the  transform  used  previously  in  Equation  2.3.  we  used 
an  -  (n  —  l/2)ir/b.  This  value  of  a„  .  however,  is  no  longer  valid  for  the  even  (in  Et ) 
modes  of  the  fin  line  we  are  calculating,  which  include  the  dominant  mode.  The  correct 
choice  is  now  a„  =  n  w/b .  This  choice  enforces  the  boundary  conditions  of  the  zero 
tangential  E  -field  and  normal  H  -field  at  Ixl  -  b.  If  we  were  interested  in  the  odd  modes  of 
the  fin  line,  the  first  definition  of  <*„.<*,  =  (n  —  l/z)yr/b  would  be  correct. 

In  general,  we  pick  a„  =  n  ir/b  when  <t>(x  )  is  even  and  Mx  )  is  odd.  where  £(x)  and 
i/»(x)  were  defined  previously  in  Equation  2.2.  This  occurs  for  even  (in  Jt )  microstrip 
modes  and  in  odd  (in  Et  )  fin  line  modes.  Furthermore,  we  pick  a,  =  (n  —  Vz)wlb  when 
d>(x)  is  odd  in  \l>(x  )  is  even.  This  occurs  in  odd  (in  7. )  microstrip  modes  and  in  even  (in 


£,  )  fin  line  modes. 
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2.7  Field  Configurations 

It  is  of  great  interest  to  plot  the  fields  due  to  the  dominant  and  higher-order  modes 
over  the  cross  section  of  the  waveguide.  This  provides  verification  that  the  boundary 
conditions  have  been  satisfied  and  offers  physical  insight  into  the  structure  of  the  modes. 
The  field  configurations  may  be  obtained  from  the  transformed  fields  obtained  in  the 
characteristic  impedance  calculation,  by  performing  the  inverse  Fourier  transform  of 
Equations  2.35.  Results  are  presented  in  the  next  section. 

2.8  Results  for  Uniform  Microstrip  and  Fin  Line 

In  this  section  we  present  numerical  results  for  the  techniques  discussed  previously 
in  this  chapter.  The  first  item  we  consider  is  the  convergence  of  the  dominant  mode 
propagation  constant  with  respect  to  the  number  of  basis  functions  and  number  of  spectral 
terms  used.  In  Table  2.1  we  show  these  calculations  and  compare  our  results  to  those  of 
Mittra  and  Itoh  [5].  From  this  table  we  can  make  a  number  of  observations.  First,  our 
values  of  0  are  in  very  good  agreement  with  those  of  Mittra  and  Itoh.  Second,  our  values 
of  0  have  converged  sufficiently  with  two  basis  functions  and  50  spectral  terms.  Note  that 
"2  basis  functions"  indicates  two  functions  for  Jz  and  two  for  Jx  .  Note  furthermore  that 
50  spectral  terms  indicates  that  all  series  were  summed  from  n  »  -49  to  50. 

Next,  we  present  a  sample  dispersion  curve  for  a  shielded  microstrip.  This  is  shown 
in  Figure  2.8.  Note  that  the  dominant  mode  is  not  cut  off.  while  the  first  higher-order  mode 
is  cut  off  below  about  20  GHz. 

Next,  we  present  data  on  the  characteristic  impedance  of  a  shielded  microstrip.  This 
is  shown  in  Figure  2.9.  The  impedance  has  been  calculated  using  both  the  V-I  and  the  V-I 
definitions,  as  discussed  in  Section  2.5,  and  the  results  are  compared  to  those  of  El-Sherbiny 
[  1 1  ].  Note  that  in  El-Sherbiny  s  paper,  the  impedance  was  calculated  for  a  shielded 
microstrip  without  side  walls.  In  order  to  account  for  this,  we  chose  in  our  calculations  to 
move  the  side  walls  far  enough  from  the  strip  to  eliminate  their  effect.  Our  results  for  the 


Table  2.1  Convergence  of  j3  with  respect  to  the  number  of  basis  functions  and  number  of 
spectral  terms.  For  these  calculations,  h  «  2  mm.  b  =  1.75  mm.  t  «  s  -  0.5  mm. 
and  er  -  9.0.  These  results  are  compared  to  those  of  Mittra  and  Itoh  [5]. 


Number  of 

Basis  Spectral 

Functions  Terms 

0(rad/m)  at 

10  GHz  20  GHz  30  GHz 

1 

25 

530.64 

1110.0 

1717.7 

1 

50 

530.50 

1109.6 

1717.1 

2 

50 

530.27 

1108.9 

1715.6 

2 

100 

530.17 

1108.7 

1715.1 

3 

250 

530.11 

1108.5 

1714.9 

Mittra  &  Itoh  [5] 

531 

1115 

1740 

V-l  definition  of  a  microstrip  agree  very  well  with  those  of  El-Sherbiny.  Although  he  does 
not  in  his  paper  specifically  identify  which  definition  of  characteristic  impedance  he  is 
using,  we  may  very  well  guess  that  it  is  the  V-I  definition. 

Next,  we  would  like  to  find  modes  in  a  fin  line.  This  is  done,  as  explained  earlier,  as 
a  check  on  the  microstrip  mode  calculations.  Results  for  the  first  three  higher-order  modes 
are  shown  in  Figure  2.10.  and  are  compared  to  those  of  Helard  et  al.  [17],  From  this  figure, 
it  is  readily  apparent  that  there  is  excellent  agreement  between  the  fin  line  modes.  This 
leads  us  to  have  a  high  degree  of  confidence  in  our  microstrip  evanescent  modes. 

Finally,  we  present  field  plots  for  the  first  five  modes  of  a  microstrip.  Since  a  given 
field  component  is  either  purely  real  or  purely  imaginary,  we  have  plotted  the  part  of  the 
field  that  is  nonzero.  In  Figures  2.10-2.13,  we  have  plotted  the  six  components  of  the 
dominant  mode  of  a  shielded  microstrip.  Note  that  these  plots  satisfy  the  boundary 
conditions  that  E-fields  tangential  to  a  conductor  and  H-fields  normal  to  a  conductor  are 
zero.  Note  also  the  correct  singularity  behavior  al  the  edges  of  the  strip.  In  Figures  2.14- 
2.17.  we  have  plotted  £,  (x  .y  )  and  £y(x  .v  )  for  the  next  four  higher-order  modes  of  the 
same  configuration.  These  are  evanescent  modes  for  the  dimensions  and  frequency  given. 
Again,  we  note  that  the  correct  boundary  conditions  and  singularities  are  observed. 


dispersion  characteristics  of  the  dominant  and  first  two  even  higher -order 
es  of  a  microstrip.  Two  basis  functions  and  100  spectral  terms  were  used. 


characteristic  impedance  of  a  microstrip.  calculated  with  two  definitions, 
basis  function  and  50  spectral  terms  were  used. 


Dispersion  curve  of  a  fin  line.  One  basis  function  and  50  spectral  terms  were 
used.  A  comparison  is  made  to  the  results  of  ilelard  et  al.  [17]. 


2.9  Conclusion 


In  this  section  we  have  demonstrated  the  capability  of  calculating  microstrip 
propagating  and  evanescent  modes.  This  extends  the  work  of  previous  authors,  as  very  few 
results  for  evanescent  modes  were  previously  available.  In  addition,  we  have  generated 
data  on  characteristic  impedance  and  fin-line  modes,  all  of  which  confirm  the  accuracy  of 
our  calculations.  Finally,  field  configurations  of  microstrip  modes  have  been  generated,  and 
these  have  been  shown  to  satisfy  the  necessary  boundary  conditions. 

All  of  the  above  suggest  a  basic  agreement  of  the  mode  calculations  with  the 
expected  results.  The  next  step  is  to  use  these  results  in  a  mode  matching  procedure  for 


CHAPTER  3 


DISCONTINUITY  CALCULATIONS 


3.1  Introduction 

The  calculation  of  discontinuities  associated  with  printed  circuits  has  been  of 
interest  for  some  time.  In  particular,  abrupt  discontinuities  in  the  widths  of  the  strips  in 
the  microstrip  and  in  the  slots  in  the  fin  line  seem  to  offer  the  best  hope  for  solution,  and 
have  been  studied  by  a  number  of  authors. 

The  earliest  attempts  at  solving  one  of  these  types  of  discontinuties  occurred  in  the 
realm  of  microstrip  discontinuities.  The  methods  used  involved  a  quasistatic  calculation  of 
a  step  discontinuity  of  strip  width  [27-30].  Although  these  methods  in  general  yielded 
good  results,  it  must  be  assumed  that  at  sufficiently  high  frequencies  the  quasistatic 
approximations  will  break  down. 

The  next  generation  of  solutions  involved  approximating  the  microstrip  as  a 
rectangular  waveguide  with  perfect  electric  conductors  on  the  top  and  bottom  walls  and 
perfect  magnetic  conductors  on  the  side  walls  [31-34],  Although  this  method  yielded  a 
frequency-dependent  solution,  it  again  is  expected  to  break  down  at  higher  frequencies  due 
to  the  approximations  of  the  model.  Furthermore,  there  are  many  printed  circuit 
discontinuities,  such  as  those  in  the  fin  line  and  strip  line,  for  which  no  waveguide  model  is 
available. 

Another  attempt  at  solving  this  type  of  problem  was  made  by  Lampe  [35],  who 
presented  a  three-region  analysis  of  strip  line  discontinuities.  Recall  that  the  strip  line  and 
fin  line  are  analogous  to  the  microstrip  because  of  the  generality  of  the  spectral  domain 
immitance  approach,  as  discussed  in  Section  2.2.  In  his  analysis,  Lampe  begins  by  breaking 
up  a  discontinuity  into  three  regions,  taking  the  boundaries  of  the  middle  region  far  enough 
away  from  the  discontinuity  on  either  side  so  that  most  of  the  evanescent  modes  have 
decayed  enough  to  become  insignificant.  Next,  he  finds  the  dominant  and  first  few  higher- 


order  modes  in  the  first  and  third  regions.  Finally,  he  formulates  an  integral  equation  that 
relates  the  magnitude  of  the  mode  coefficients  in  the  first  and  third  regions  to  the  unknown 
current  on  the  strip  in  the  second  region.  The  solution  to  this  integral  equation  yields  the 
desired  reflection  and  transmission  coefficients.  This  method  has  limitations  also,  because  it 
is  formulated  for  an  isotropic  medium.  The  strip  is  assumed  to  be  printed  on  a  thin 
dielectric  substrate;  so  thin  that  it  is  felt  to  have  little  effect.  Thus,  the  entire  problem  is 
formulated  in  a  homogeneous  region  whose  dielectric  constant  is  the  "volume  average" 
dielectric  constant.  This  turns  out  to  be  a  small  perturbation  from  the  free-space  dielectric 
constant  of  one.  This  method  is  appropriate  for  low  frequencies  in  the  strip  line,  but  it 
seems  likely  to  break  down  for  higher  frequencies.  In  addition,  it  is  not  appropriate  for  a 
microstrip  configuration,  in  which  the  dielectric  material  plays  an  important  role.  Finally, 
this  method  can  not  easily  be  modified  to  account  for  the  dielectric,  since  the  Green's 
function  in  the  second  region  would  then  become  too  complicated. 

The  most  recent  attempts  at  solving  this  type  of  problem  have  occurred  in  the  area 
of  fin  line  discontinuities.  The  methods  used  in  these  cases  all  involved  finding  a  number  of 
modes  for  each  guide  by  using  a  spectral  Galerkin  method.  These  modes  were  then  used  in 
either  a  mode  matching  procedure  [  17.36.37]  cr  in  an  iteration  procedure  [38]  to  equate  the 
fields  tangent  to  the  plane  of  the  discontinuity. 

The  results  generated  in  these  mode  matching  papers  were  not  conclusive,  but  they 
did  suggest  that  further  study  would  be  in  order,  and  that  the  method  might  be  applicable 
to  a  shielded  microstrip  as  well  as  to  the  fin  line.  This  will  be  the  topic  of  study  for  this 
chapter. 

3.2  Mode  Matching 

The  structure  to  be  analyzed  is  shown  in  Figure  3.1a.  It  involves  an  abrupt 
discontinuity  in  the  strip  width  of  the  center  conductor  of  shielded  microstrip,  whose  cross 
section  is  shown  in  Figure  2.1. 
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We  begin  the  analysis  by  generating  the  propagation  constants  and  field 
configurations  for  the  dominant  and  first  few  higher-order  modes  of  each  of  the  two 
waveguides.  This  was  discussed  in  Chapter  2.  Next,  we  express  the  fields  in  waveguide  A 
as  a  sum  of  the  waveguide  modes 

~Ea  (x  ,y  j.  )  =  £  a;  eai  (x  .y  )  e  ~>J‘ 2  (3.1a) 

i=l 

Ha  (x  .y  .2  )  =  £  a,  Rai  (x  .y)e  *  (3.1b) 

i  =i 

where  the  a,  are  mode  coefficients.  yai  are  the  complex  propagation  constants  of  the  ith 
mode,  and  fa,  and  Rai  are  the  vector  mode  f unctions  of  the  ith  mode,  where 

ea,  (x.y)-  eatx  ( x  ,y  )x  +  ea,y  (x  .y  )y  (3.2a) 

K<,i  (x  .y  )  =  h (x  ,y  )x  +  h„ix  (x  .y  )y  (3.2b) 

Similar  expressions  are  formulated  for  waveguide  B  by  replacing  ’a’  with  *b"  in  the  above 
equations.  If  we  assume  the  discontinuity  occurs  at  z  -  0.  we  may  equate  the  field 
components  tangential  to  the  interface  at  z  •  0  as 

oo  ao 

(1  +p)ea  lU  •>•)+£  a,  «a,  (x.y)  -  Z  *>>  U  .y  )  (3.3a) 

i  ~2  i*l 

( l-p)/ra  i(x  .y  )  -  £  a,  Rai  (x  .y  )  =  £  b,  Rb,  (x  .y  )  (3.3b) 

i*2  i=l 

This  is  the  equation  that  is  to  be  solved  for  the  mode  coefficients  a,  and  b,  and  for  the 
reflection  coefficient  p. 

To  solve  this  equation,  we  take  the  inner  products  of  Equation  3  3a  with  ^.(x  .y  ), 
and  the  inner  product  of  Equation  3.3b  with  (x.y).  The  inner  product  is  defined  as 

h,b,  -  f  f  *a,(*  -y)  *  Kb,  (x  .y  )  z  dx  dy  (3.4) 

s 

It  is  calculated  most  simply  in  the  spectral  domain,  in  a  manner  analogous  to  the  power 


calculation  described  in  Section  2.5.  Thus,  the  integration  over  x  becomes  a  sum  over 
spectral  terms,  and  the  integration  over  y  is  analytic  in  the  spectral  domain.  After  taking 
these  inner  products,  we  find 


P  1  am  X  ^ atom  ^  ^ btam  lam  ^  —1......Z,  (3. 5a) 

i  =2  i*l 

L  L 

P  Ibma  1  ^  f bmai  22  ?  bmbi  I  m  —  1 . Z.  (3.5b) 

i=2  i=l 

This  equation  can  also  be  written  in  the  form 

Ax  =  b  (3.6) 


where  A  is  a  2L  X  2L  matrix  of  inner  products,  x  is  the  2L  X  1  vector  containing  the 
unknown  mode  coefficients,  and  b  is  a  known  2L  X  1  matrix  of  inner  products.  This  can  be 
solved  by  a  standard  linear  equation  solution  routine. 

The  solution  to  this  matrix  equation  contains  within  it  the  reflection  coefficient,  p, 
and  the  magnitude  of  the  transmitted  wave,  b  t.  From  p.  we  may  calculate  an  equivalent 
circuit  of  the  form  shown  in  Figure  3.2.  If  we  let 

P  =  Pr  +  j  P,  (3.7) 


where  pr  and  p,  are  both  real,  then 


Zy 


Z2 


Pr  +  1 

2  +  p,2 

i-p; 

-p,2 

(3.8a) 


Y  y  =  ojCZ  ,  = 


-2p, 

Pr  +  1  2  +  Pi2 


(3.8b) 


These  parameters  are  useful  because  they  can  be  checked  with  other,  more  approximate 
methods.  Therefore,  we  expect  Z.v  *=  1*  /Z„  where  and  Zx  are  the  characteristic 
impedances  of  the  second  and  first  waveguides,  as  calculated  in  Section  2.5.  Furthermore, 
we  may  compare  our  values  of  T.v  to  those  generated  by  quasistatic  analysis  [29].  We 


should  expect  only  approximate  agreement,  since  we  are  now  using  a  frequency  dependent 
method  to  generate  the  solution. 

A  second  check  we  can  make  on  our  final  answer  is  to  verify  that  all  real  power  is 
accounted  for.  This  is  checked  by  showing  that 

IS„I2  +  IS21I2  =  1  (3.9) 

While  this  condition  is  not  sufficient  to  guarantee  a  good  solution,  it  is  necessary,  and 
should  be  verified. 

This  concludes  the  setup  of  the  mode  matching  procedure.  Before  we  present 
calculations,  there  are  a  number  of  related  issues  that  need  to  be  dealt  with.  These  are 
discussed  in  the  following  three  sections. 

3.3  Orthogonality  of  Inner  Products 

Upon  examination  of  Equation  3.5  we  observe  a  number  of  inner  products  of  the 
form  IaiaJ  where  i  j.  We  expect  the  standard  orthogonality  relationship  to  hold  for  the 
normalized  modes  of  uniform  waveguides  with  perfectly  conducting  side  walls  [39] 

Alia;  =  o  i  9s  )  (3.10) 

Thus,  we  are  left  with  the  problem  of  deciding  whether  or  not  to  retain  in  the  matrix 
equation  the  inner  products  that  are  theoretically  equal  to  zero.  It  must  be  kept  in  mind 
that  our  modes  have  been  calculated  only  to  a  finite  accuracy,  dictated  by  the  numerical 
accuracy  of  our  methods.  If  the  accuracy  in  the  mode  functions  is  good,  then  retaining  the 
cross  terms  will  not  matter.  If.  on  the  other  hand,  there  is  a  small  error  in  the  mode 
functions,  then  it  would  seem  preferable  to  retain  the  cross  terms  rather  than  discard  what 
may  be  useful  numerical  information.  This  is  the  approach  adopted  for  these  calculations. 
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3.4  Condition  Number  of  the  Matrix 

In  formulating  the  mode  matching  solution,  we  obtain  a  matrix  equation  of  the 
form  Ax  -  b.  A  useful  parameter  associated  with  this  equation  is  the  condition  number  of 
the  matrix  A.  It  turns  out  that  if  the  matrix  A  has  a  large  condition  number,  it  is  very 
difficult  to  obtain  accurate  solutions  for  the  unknown  column  vector  x. 

There  are  a  number  of  methods  that  may  be  used  to  calculate  the  condition  number 
of  a  matrix,  depending  on  the  norm  of  the  matrix  chosen  for  the  problem  [40.41],  A 
reasonable  definition  of  condition  number  is 


cond  (A  )  = 


l\ 

U 


mu  ^ 
mm  ^ 


(3.11) 


where  the  \'s  are  the  maximum  and  minimum  eigenvalues  of  the  matrix  AH  A  .  and  AH  is 
the  transposed  complex  conjugate  of  A.  If  this  condition  number  is  large,  we  have  difficulty 
in  solving  the  matrix  equation,  because  small  errors  in  the  matrix  elements  generate  large 
errors  in  the  solution  for  the  unknown  vector  x.  This  definition  of  condition  number  is 
used  in  the  calculations  that  follow,  in  order  to  check  the  stability  of  the  matrix  equation 
solution. 


3-5  Matrix  Theory  for  Cascaded  Discontinuities 

In  the  case  where  there  are  multiple  discontinuities,  we  need  a  method  of  keeping 
track  of  a  number  of  modes  between  one  discontinuity  and  the  next.  We  consider  two 
cases.  First,  we  study  a  symmetrical  double  step,  as  shown  in  Figure.  3.1b.  This  is  a 
special  case  of  the  second  case  we  will  study,  that  of  N  discontinuities  each  spaced  an 
arbitrary  distance  from  the  previous  discontinuity.  An  example  of  this  is  shown  in  Figure 
3.1c. 

In  the  case  of  an  asymmetrical  double  step,  we  may  take  advantage  of  svmmetrv 
properties  to  greatly  simplify  the  problem.  Thus,  instead  of  launching  the  dominant  mode 


n 


from  the  left  with  a  magnitude  of  one,  we  add  the  results  of  launching  modes  of  equal 
magnitude  and  either  equal  or  opposite  phase,  as  shown  in  Figure  3.3.  These  cases  are 
equivalent  to  placing  magnetic  and  electric  conductors  at  the  center  of  the  waveguide  B. 
The  solution  for  each  of  these  cases  may  easily  be  adapted  from  the  solution  for  the  single 
discontinuity,  since  the  waveguide  modes  are  reflected  from  the  electric  or  magnetic 
conductor  without  coupling  to  other  modes.  We  adapt  results  from  [36]  that  result  in  a 
slight  modification  of  Equation  3.3.  Thus,  at  z  -  0 

( 1/2  +  pVa  i(*  .y  )  +  £  a,  eai  (x  .y  )  =  £  b>  ( 1  +  Sbii  )eb,  (x  ,y  )  (3.12a) 

i  =2  i  *1 


( 1/2  -  p)/?a  i(x  .y  )  -  £  a,  Kai  (x  ,y  )  =  £  b,  ( 1  -  Sbil  )Kb,  (x  .y  ) 

i  =2  i  =1 


(3.12b) 


where 


Sbu  — 


-2  >»,' 


e  for  electric  conductor  at  z  =  / 


-e  2y>’'  for  magnetic  conductor  at  z  ~  l 


(3.13) 


and  yb,  is  the  complex  propagation  constant  of  the  ,th  mode  in  waveguide  B.  and  21  is  the 
distance  between  the  discontinuities.  The  above  equation  may  be  solved  in  a  manner 
analogous  to  Equation  3.3  by  taking  appropriate  inner  products,  as  described  in  Section  3.2. 

If.  on  the  other  hand,  we  do  not  have  the  luxury  of  having  favorable  symmetry 
properties,  we  have  to  solve  the  more  general  case  of  N  discontinuities,  each  separated  by  a 
length  of  transmission  line.  In  this  case  we  will  have  to  cascade  the  generalized  S- 
parameters  of  each  of  the  discontinuties.  and  the  transmission  lines  that  separate  them,  to 
form  an  S  matrix  for  the  entire  structure.  The  method  that  follows  is  adapted  from  a 
method  by  Hall  [42]. 

Let  us  begin  with  a  definition  of  generalized  S-  and  T-parameters.  For  the  arbitrary 


circuits  shown  in  Figure  3.4.  we  define  generalized  T-  and  S-parameters  as 


Electric 

Conductor 


1/2 


Magnetic 

Conductor 


Symmetrical  discontinuities,  and  the  method  ol  taking  advantage  ol'  the 
symmetry  by  symmetrical  and  anlisy  m metrical  excitation  ol  the 
discontinuity. 
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*2 

T  i, 

T  12 

“1 

a2 

T  21 

T  22 

*1 

b  i 

^11 

■$12 

“l 

b  2 

$21 

S22 

“2 

(3.14a) 


(3.14b) 


In  these  equations,  a,  and  b ,  represent  L  X  1  column  vectors  representing  mode  coefficients 
for  the  L  modes  either  incident  or  reflected  from  the  circuit.  In  addition,  terms  such  as  Su 
represent  an  L  X  L  generalized  scattering  matrix  for  L  modes  in  each  guide. 

We  generate  generalized  S-matrices  for  discontinuities  in  the  manner  of  Section  3.2. 
but  keeping  in  mind  that  we  now  have  to  run  cases  for  all  L  possible  modes  incident  from 
the  left  at  each  discontinuity.  Care  must  be  taken  with  the  normalization  of  these  S- 
parameters,  since  the  modes  generated  were  not  normalized.  The  correct  normalization 
gives,  for  example 


S21(i  ,j ) 


a, 


f  febJ  X  Kbj  d  S 
f  f  ?oi  d~S 


(3.15) 


Once  the  generalized  scattering  parameters  are  normalized,  they  must  be  converted  to  T- 
parameters.  The  conversions  are 


*^21 

~  $  22^  ll'^  11 

$22$  12* 

-Sw'Sn 

Sf2* 

-Ti2lT2i 

T22 

Tn 

~  T  \2T  22  T  21 

Tx2Tf2' 

(3.16a) 


(3.16b) 


Next  we  generate  the  T-matrices  of  individual  transmission  lines.  These  are 


Tuu.n 


•  =j 

else 


(3.17a) 
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parameters,  Z.v  and  r  ,v  •  are  shown  in  Table  3.2  for  calculations  with  up  to  six  modes  in 
each  waveguide.  We  note  here  that  Z.v  converges  quickly  to  the  values  expected  from 
transmission  line  theory.  The  junction  capacitance,  however,  never  seems  to  converge  with 
the  small  number  of  modes  we  have  used.  We  are  limited  somewhat  in  the  number  of 
modes  we  can  use  because  the  cost  of  the  calculation  begins  to  become  quite  large  with  more 
modes,  and  because  the  spectral  Galerkin  technique  tends  to  break  down  for  large 
imaginary  propagation  constants. 

Next,  we  perform  a  number  of  checks  on  these  calculations.  These  include  a  power 
check  and  a  calculation  of  the  condition  number  of  the  A  matrix.  These  data  are  shown  in 
Table  3.3.  From  this  table  we  verify  that  all  the  power  is  accounted  for  to  a  reasonable 
numerical  accuracy.  This,  as  stated  earlier,  is  a  necessary,  but  not  sufficient  condition  to 
guarantee  the  accuracy  of  a  solution.  The  condition  number  of  the  matrix  may  be  more 
indicative  of  what  is  going  on.  The  condition  numbers,  which  are  around  750.  suggest  that 
the  matrix  is  ill-conditioned.  As  a  general  guideline,  we  consider  a  condition  number  greater 
than  about  100  to  indicate  a  problem  in  the  condition  number  of  the  matrix.  Thus,  unless 
the  matrix  elements  are  computed  very  accurately,  we  have  to  expect  a  problem  in 
generating  highly  accurate  solutions  to  the  matrix  equation. 

The  issue  of  condition  number  turns  out  to  be  of  a  very  central  importance  in  these 


Table  3.2.  Values  of  the  circuit  elements  in  Figure  3.2  for  the  dimensions  in  Table  3.1. 

These  calculations  were  made  with  5  basis  functions  and  200  spectral  terms. 


#  of  Modes 

Pr 

P, 

Y.y 

2,v 

3 

-.269978 

-.0014 

.0060 

.575 

4 

-.270028 

-.0011 

.0043 

.575 

5 

-.270368 

-.0038 

.0146 

.574 

6 

-.280095 

-.0060 

.0233 

.562 

Expected  from 
Quasi  Static  0"'v ) 

.0172 

.684 

and  TEM  (Z.v ) 

Approaches 

Table  3.3.  Power  check  and  condition  number  for  the  discontinuity  whose  dimensions 
were  given  in  Table  3.1. 


Number 
of  Modes 

IS  ii  1 2  +  IS2,I2 

Condition 

Number 

3 

1.0016 

754 

4 

1.0016 

759 

5 

1.0006 

752 

calculations.  Since  the  imaginary  part  of  the  reflection  coefficient  is  several  orders  of 
magnitude  less  than  the  real  part,  the  reflection  coefficient  as  a  whole  must  be  calculated  to 
a  high  degree  of  precision  in  order  to  get  a  meaningful  p, .  and  hence  YN.  When  the 
condition  number  is  high,  it  is  very  difficult  to  obtain  accurate  results  without  a  very  high 
degree  of  numerical  precision  in  the  matrix  elements. 

Let  us  now  consider  the  accuracy  of  the  matrix  elements.  The  simplest  way  of 
checking  this  is  to  consider  the  degree  to  which  the  orthonormality  conditions  of  the 
waveguide  modes,  as  expressed  in  Equation  3.10.  were  satisfied.  If  the  orthogonality 
conditions  are  satisfied  well,  then  we  may  concede  the  possibility  of  an  accurate  solution  to 
the  discontinuity  problem  despite  the  ill-condition  of  the  matrix.  These  inner  products 
appear  in  Table  3.4  for  the  first  five  modes  of  waveguides  A  and  B.  From  this  table,  we  see 
that  the  cross  terms  are  well  behaved  for  the  lower-order  terms,  but  the  higher-order  cross 
terms  tend  to  become  increasingly  large.  Thus,  it  is  difficult  to  claim,  based  on  these  inner 
products,  that  the  matrix  elements  are  accurate  enough  to  overcome  the  large  condition 
number  of  the  A  matrix. 

Finally,  we  present  data  on  the  number  of  basis  functions  and  spectral  terms  that 
are  required  to  give  accurate  propagation  constants.  This  may  be  used  to  answer  a  possible 
objection  that  an  insufficient  number  of  basis  functions  and  spectral  terms  was  used  in  the 
discontinuity  calculations.  The  data,  shown  in  Table  3.5.  indicate  that  the  propagation 
constant  of  the  fifth  mode  in  waveguide  A  is  sufficiently  converged  with  only  two  basis 
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Table  3.4.  Inner  product  calculations  for  the  first  five  modes  of  waveguides  A  and  B, 
whose  dimensions  were  given  in  Table  3.1.  These  were  calculated  with  5  basis 
functions  and  200  spectral  terms. 


i  1 

1.000 

0.000 

0.000 

0.035 

0.013 

0.000 

1.000 

0.000 

0.009 

0.003 

0.000 

0.000 

1.000 

0.048 

0.018 

0.000 

0.000 

0.000 

1.000 

0.010 

0.000 

0.000 

0.001 

0.521 

1.000 

1/6*;  1 

1.000 

0.003 

0.046 

0.001 

0.080 

0.000 

1.000 

0.053 

0.000 

0.009 

0.000 

0.000 

1.000 

0.001 

0.086 

0.000 

0.000 

0.113 

1.000 

0.031 

0.000 

0.000 

0.031 

0.000 

1.000 

functions  and  50  spectral  terms.  Thus,  since  we  used  five  basis  functions,  and  200  spectral 
terms  in  our  discontinuity  calculations,  it  seems  likely  that  we  have  used  a  sufficient 
number  of  each. 

We  conclude,  therefore,  that  the  mode  matching  technique  is  useful  only  in 
obtaining  good  approximations  to  the  circuit  parameters.  It  seems  unlikely  that  the 
accuracy  of  this  method  can  be  forced  to  the  point  where  an  accurate  junction  capacitance 
can  be  calculated,  for  a  number  of  reasons.  First,  the  imaginary  part  of  the  reflection 
coefficient  is  very  small  compared  to  the  real  part,  and  it  is  difficult  to  calculate  a  small 


Table  3.5.  Variation  of  a  propagation  constant  with  the  number  of  basis  functions  and 
spectral  terms.  The  mode  calculated  is  the  fifth  mode  of  waveguide  A.  whose 
dimensions  were  given  in  Table  3.1. 


Number  of 

Basis  Spectral 

Functions  Terms 

s(rad/m) 

1 

25 

14901.0 

2 

50 

14897.1 

2 

100 

14894.5 

3 

50 

14897.6 

3 

100 

14894.7 

5 

250 

14893.3 

I 
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quantity  in  the  shadow  of  a  larger  effect.  Second,  the  condition  number  of  the  matrix 
indicates  an  instability  in  the  matrix  that  can  only  be  overcome  if  the  matrix  elements  are 
calculated  to  a  high  degree  of  accuracy.  Finally,  the  inner  product  calculations  suggest  that 
the  matrix  elements  can  only  be  calculated  to  a  finite  degree  of  accuracy  and  that  the 
spectral  Galerkin  technique  can  not  be  pushed  beyond  this  point  for  modes  of  high  order. 

With  these  thoughts  in  mind,  let  us  now  turn  to  other  discontinuities  to  calculate. 
Although  we  have  not  achieved  a  high  degree  of  accuracy  in  the  calculations  for  the  single 
discontinuity,  we  have  demonstrated  that  the  method  generates  a  good  approximation  for 
the  solution.  Thus,  there  may  well  be  reason  to  consider  other  types  of  discontinuities,  and 
results  for  these  are  presented  in  the  section  that  follows. 

3.7  Results  for  Other  Discontinuities 

The  next  configuration  we  would  like  to  study  is  the  symmetrical  double  step 
discontinuity.  This  is  shown  in  Figure  3.5,  and  the  theory  was  presented  earlier  in  Section 
3.5.  Results  are  presented  for  a  typical  case  at  two  different  frequencies  in  Tables  3.6  and 
3.7.  These  calculations  were  made  with  up  to  four  waveguide  mode  functions  in  each 
waveguide.  Upon  examination  of  these  results,  we  find  the  reflection  coefficient  has 
converged  nicely  within  four  modes  to  a  result  that  is  similar  to  that  expected  from  the 
transmission  line  theory.  We  note,  furthermore,  that  the  results  for  one  mode  are  similar 
to  that  for  four  modes,  so  in  the  future  we  need  to  use  only  a  single  mode  for  our 
calculations. 

In  the  next  two  tables.  Tables  3.8  and  3.9,  we  present  data  for  several  of  these  cases 
over  a  range  of  frequencies  and  for  various  values  of  sb .  We  compare  them  to  experimental 
data,  which  was  generated  by  U.  Feldman  [43],  and  to  results  generated  by  the  transmission 
line  theory.  Based  on  the  data  in  these  tables,  we  observe  that  the  data  calculated  by  the 
mode  matching  technique  provide  a  slightly  belter  fit  to  the  experimental  data  than  the 
results  generated  by  the  transmisison  line  theory. 


Table  3.6. 


Propagation  constants  for  the  first  through  fourth  modes  of  a  symmetrical 
double  step  discontinuity  in  strip  width,  and  results  for  Sa-  For  these  calcu¬ 
lations.  h«5.08  mm,  b-6.096  mm.  t=0.7874  mm.  er  -2.2.  Z-1.0  cm.  and 
freq-8.010  GHz. 


^(rad./mm) 

02 

03 

04 


«  1.17  mm 
,  -  51.85  n 
230.056 
-jl87.051 
-j647.067 
-j727.228 


sb  ■  2.41  mm 
Zpt,  -  3.60  n 
234.723 
-jl86.281 
-j640.983 
-j749.948 


Table  3.7. 


Propagation  constants  for  up  to  the  fourth  mode  of  a  symmetrical  double  step 
discontinuity  in  strip  width,  and  results  for  Su.  The  configuration  is  that  of 
Table  3.6.  except  that  the  freq  *  12.02  GHz. 


341.625 

355.340 

55.8937 

55.8765 

-j701 .537 

-j610.145 

-j723.682 

-j904.769 

Number 
of  Modes 

S„  ( dB  ideg  ) 

1 

-8.90  l  -140.6 

2 

-8.90  i  -140.6 

3 

-8.46  l  -140.3 

4 

From  Transmission 

-8.45  l.  -140.3 

Line  Theory  Expect 

-10.8  l  -139.4 

Table  3.8. 


Dominant  mode  propagation  constant  and  characteristic  impedance  as  a  func¬ 
tion  of  frequency  and  center  strip  half  width  (s).  For  these  calculations,  h  « 
5.08  mm,  b  -  6.096  mm.  t  -  0.7874  mm,  and  €r  -  2.2. 

/3  (rad/m)  and 
Z0  (H)  for 


FREO  (GHz) 

s  -  2.41 

s  *  1.65 

s-  1.17 

s  -  .026 

s  -  .353 

.990 

28.6661 

34.52 

28.43 

42.26 

28.1895 

51.74 

27.7751 

73.67 

27.4819 

97.03 

2.025 

58.6799 

34.50 

58.1802 

42.23 

57.69 

51.72 

56.8345 

73.65 

56.22% 

97.01 

4.005 

116.374 

34.45 

115.330 

42.19 

114.319 

51.68 

112.563 

73.61 

111.330 

%.97 

8.010 

234.723 

34.60 

232.338 

42.35 

230.056 

51.85 

226.170 

73.81 

223.470 

97.19 

12.015 

355.340 

35.21 

351.391 

43.01 

347.625 

52.56 

341.250 

74.60 

336.825 

98.06 

Table  3.9. 


Su  for  the  symmetrical  double  step  discontinuity  shown  in  Figure  3.4.  as  a 
function  of  sb  and  frequency.  For  these  calculations.  sa  -  1.17  mm.  I  -  1.0 
cm.  and  all  other  dimensons  are  as  in  Table  3.8.  The  first  number  for  each  case 
is  calculated  by  mode  matching,  the  second  by  the  transmission  line  theory, 
and  the  third  number  is  from  experimental  results  of  Feldman  [43]. 

_ _ 5  n  (dB  l  deg)  for 


FREO  (GHz) 

sh  •  2.41  mm 

sh  «  1.65  mm 

sh  -  .626  mm 

sh  «  .353  mm 

.990 

-10.8/.  -126 
-13.2/. -125 
-11.2 

-17.6/. -123 
-19.2/. -123 
-18.3 

-13.9/.  56 
-14.6  4  57 
-13.7 

-9.1  4  53 
-9.6  4  54 
-9.4 

2.025 

-6.8  L  -160 
-8.9/. -159 
-6.9 

-13.1  l  -157 
-14.7/. -157 
-13.3 

-9.5  L  23 
-10.1  4  24 
-9.5 

-5.2  4  21 
-5.7  4  22 
-5.4 

4.005 

-8.6  L  140 
10.8  L  139 
-8.9 

-14.9  L  138 
-16.5  L  138 
-15.3 

-10.7  4  -37 
-  -11.4 /.-37 
-11.1 

-6.1  4-32 
-6.6  4  -33 
-6.6 

8.010 

-6.4/. -179 
-8.3/. -179 
-6.4 

-12.5  l  -176 
-14.0/. -176 
-12.6 

-8.9  L  10 
-9.5  4  10 
-9.4 

-4.9  4  11 
-5.3  4  12 
-5.2 

12.015 

-8.9  l  -141 
-10.8/. -139 
-8.8 

-15.9/. -133 
-17.4/.  -133 
-16.4 

-14.3/.  57 
-14.8/.  57 
-14.7 

-10.7  4  59 
-11.1  46O 
-10.3 

Next,  we  consider  a  nonsymmetrical  double  discontinuity,  shown  in  Figure  3.6. 
This  configuration  may  be  used  to  simulate  a  linear  taper.  The  results  for  this  structure 
were  generated  with  the  matrix  theory  given  in  Section  3.5.  and  are  given  in  Table  3.10.  In 
this  table  are  the  propagation  constants  and  characteristic  impedances  of  the  three  lines,  and 
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the  S-parameters  of  the  discontinuity  referenced  to  planes  Nos.  1  and  2  as  shown  in  Figure 
3.6.  In  addition.  Table  3.10  has  a  comparison  to  transmission  line  theory  and  a  power  check 
of  the  mode  matching  results.  The  results  indicate  an  overall  agreement  with  the 
transmission  line  approach,  although  it  is  difficult  to  say  which  approach  is  more  accurate. 

Finally,  we  consider  a  single  discontinuity  in  the  dielectric  constant,  strip  width, 
and  substrate  thickness,  a  diagram  of  which  is  shown  in  Figure  3.7.  This  case  is  one  that 
might  be  expected  to  occur  when  a  microstrip  printed  on  gallium  arsenide  must  mate  with  a 
microstrip  printed  on  duroid  board.  Because  of  the  difference  in  dielectric  constants,  it  will 
be  necessary  to  have  different  line  widths  to  maintain  a  50  ohm  line  in  each  section.  Some 
typical  S-parameter  data  for  this  configuraton  are  shown  in  Table  3.11.  along  with  data  for 
a  power  check  and  the  condition  number  of  the  matrix.  These  data  are  difficult  to  interpret 
since  the  condition  numbers  are  very  large,  and  since  the  power  check  is  off  by  about  0.06. 
Both  transmission  lines  are  50  ohm  lines,  so  we  expect  a  reflection  coefficient  of  zero  from 
simple  transmission  line  analysis,  while  we  calculate  a  reflection  coefficient  of  about  -10  dB. 
This  -10  dB  reflection  corresponds  to  a  value  of  0.1  in  the  power  check.  Since  the  power 
check  is  off  by  0.06.  it  is  difficult  to  get  a  feel  for  the  accuracy  of  these  results.  If  we 
ignore,  however,  the  higher -order  modes,  and  look  only  at  the  results  when  a  single  mode  is 
used,  we  still  have  a  reflection  coefficient  of  about  -10  dB.  but  now  the  power  check  and 
condition  number  of  the  matrix  are  both  satisfactory.  It  appears,  therefore,  that 
experimental  work  will  be  required  to  verify  these  calcinations. 

3.8  Conclusion 

In  this  chapter  we  have  analyzed  a  number  of  discontinuities  by  using  a 
combination  of  the  spectral  Galerkin  technique  to  generate  modes  and  mode  matching  to 
find  the  scattering  parameters  of  the  discontinuity.  In  general,  our  results  have  been  close 
to  what  was  expected,  but  it  proves  difficult  to  use  this  technique  to  give  highly  accurate 
results.  The  factors  that  limit  the  accuracy  include  the  small  number  of  waveguide  modes 
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Table  3.10.  Results  for  the  taper  simulation  shown  in  Figure  3.6.  For  these  calculations, 
h  «  5.08  mm.  b  -  6.096  mm.  t  -  0.7874  mm.  er  -  2.2,  and  l  -  1.0  cm. 


' 1 

10  GHz 

12.5  GHz 

15.0  GHz 

sa  -  1.17  mm 

(rad/m) 

EISlHi 

435.9 

Zaz  (fl) 

■ 

53.28 

sb  «  1.65  mm 

291.2 

Kf^EfWIiiW 

440.8 

Z-ob 

42.63 

43.67 

r0  =  2.41  mm 

Pc 

294.36 

269.9 

445.9 

Zx 

34.86 

35.30 

35.79 

Transmission 

S ii 

-14.95  /  -154 

-19.76/.  122 

-15.81  /  -146 

Line 

21 

-0.14/.  26 

-0.05  /  -59 

-0.12  / -146 

Theory 
dB  /  deg 

$22 

-14.95  /  26 

-19.76/  -59 

-15.81  /  34 

Mode 

5  11 

-12.78  / -152 

-17.61  /  114 

-12.94/ -143 

Matching 

•$21 

-0.23  /  26 

-0.76  /  -58 

-0.18  / -146 

s  22 

-12.78  /  23 

-17.61  / -51 

-13.94  /  31 

IS,,I2  +  IS21I2 

1.00000 

0.999998 

0.999999 

Table  3.11.  Results  for  a  step  discontinuity  in  €r  .  t.  and  s.  For  these  calculations  ha  = 
0.889  mm.  hb  -  0.953  mm.  b  »  1.27  mm.  ta  «  0.127  mm.  tb  «•  0.191  mm.  sa  » 
0.042  mm.  sb  =  0.298  mm.  =*  12.3.  erb  =  2.2.  freq.  «  20  GHz.  These  calcu¬ 
lations  were  made  with  2  basis  functions  and  50  spectral  terms,  and 
ZM  =  =50.0  n. 


Number 
of  Modes 

Su  ( dB  /  deg  ) 

S2j  (dB  /  deg  ) 

ISji  1 2  +  IS21I2 

Condition 

Number 

1 

-9.92  /  180 

-.46/0.0 

1.00000 

26 

2 

-11.69/ -175 

-.08  /  -2.2 

1.066 

6240 

3 

-1 1.67  /  -1 75 

-.01  /  -2.0 

1.066 

6790 

4 

-1 1.67  /  -172 

-.02  /  -3.2 

1.063 

6800 

61 


one  can  calculate  with  the  spectral  Galerkin  technique,  the  ill-condition  of  the  set  of  matrix 
equations,  and  the  necessity  of  calculating  a  small  junction  capacitance  in  the  shadow  of  a 
comparatively  large  reflection  coefficient. 

If  this  method  is  to  be  refined  to  yield  more  accurate  results,  we  have  to  find  a  way 
to  generate  a  large  number  of  very  accurate  waveguide  modes.  In  the  next  chapter,  we 
consider  an  alternate  method  of  generating  these  modes,  in  an  effort  to  increase  the  accuracy 
of  the  mode  matching  method.  Let  us  now  turn  our  attention  to  the  singular  integral 
equation  technique. 
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CHAPTER  4 

THE  SINGULAR  INTEGRAL  EQUATION  TECHNIQUE 

4.1  Introduction 

The  singular  integral  equation  technique  is  an  alternate  method  of  finding  modes  in 
printed  circuit  waveguides.  This  method  was  first  suggested  by  Mittra  and  Itoh  [5],  and  has 
since  been  expanded  upon  by  Omar  and  Schunemann  [44.45].  Its  chief  advantage  over  the 
spectral  Galerkin  method  is  its  very  high  degree  of  numerical  efficiency,  requiring,  by 
comparison,  very  little  computer  time.  Its  disadvantage,  however,  is  that  it  is  a  much  more 
difficult  method  to  derive  and  implement,  requiring  a  great  deal  of  effort  to  derive  the 
matrix  equation. 

In  this  chapter,  we  outline  the  work  of  Mittra  and  Itoh  and  expand  upon  it  by 
calculating  a  larger  number  of  terms  in  certain  series  as  well  as  a  larger  matrix  size.  In 
addition,  we  compare  these  results  to  those  obtained  with  the  spectral  Galerkin  technique 
presented  earlier. 

4.2  Overview  0f  the  Method 

In  this  section  we  outline  the  singular  integral  equation  technique  as  applied  to  the 
microstrip,  and  as  explained  in  [5].  Since  the  method  is  somewhat  involved,  we  omit  a 
number  of  details  here  that  are  included  in  [5],  In  spite  of  this,  we  include  sufficient  detail 
to  demonstrate  the  reasons  for  the  increased  efficiency  of  the  method. 

Consider  the  shielded  microstrip  structure  shown  in  Figure  4.1.  We  may  write  the 
fields  in  the  two  regions,  denoted  by  i.  as  a  sum  of  TE  and  TM  components.  This  leads  to 

k  2  -  Q2 

E,,  -  j-L— - (4.1a) 

E,1:1  *  V,  <£,u’  e~‘ 


» 


(4.1b) 


#,<'>=  '  £  x 

P 


(4.1c) 


for  the  TM  fields  and 


k  2  -  Q2 

Hzi  =  j 


(4.2a) 


Et\h  >  =  -  £  XV,  il'i1'1  ’e-' *2 


(4.2b) 


H,\h)  =  V,  0/A)e“^2 


(4.2c) 


for  the  TE  fields.  In  these  equations  /3  is  the  propagation  constant  in  the  £  direction.  <Jr/f * 
and  <Jr/*  *  are  scalar  potentials  for  the  E  and  H  fields,  respectively,  and  the  subscript  t 
indicates  the  transverse  direction.  The  scalar  potentials  each  satisfy  a  two-dimensional 
Helmholtz  equation  as  in  Equation  2.2  of  the  spectral  Galerkin  approach,  where 
k,  =  (oy/fx0  €,"  is  the  propagation  constant  of  the  medium  in  regions  1  and  2.  By  satisfying 
the  boundary  conditions  on  the  side  and  top  walls,  we  expand  the  scalar  potentials  as 


'hi' )  =  X  Anr )  sinh  a„(,,y  cos  k„  x 


(4.3a) 


hi' )  =  £  Bn{' )  sinh  a„(2)  {h  —  y)  cos  kn  x 


(4.3b) 


hih }  =  £  An(h  ;  cosh  a„a)y  sin  icn  x 


(4.3c) 


hih  >  =  £  B ’„(A  >  cosh  aj2)  (h  -  y  )  sin  kn  x 


(4.3d) 


where 


(n  —  V2)V 


2 


and  k0  —  fi0  €0  . 


We  may  now  apply  the  boundary  conditions  at  y  -  d,  which  is  the  interface  of  the 


two  regions,  to  obtain 


Z  An(t  >  cos  kn  x  -0. 

'  n  «1 


Z  Af'kn  PAfi)cosknx 


0  <  x  <  t 


(4.6a) 


-  Z  Anh  )  K  T„  (0)  cos  kn  X  =  0.  t  <  X  <  L  (4.6b) 


£  Anu  1  kn  sin  kn  x 


~  Z  Anh  )  K  sin  k„  x  =  0.  0  <  x  <  t 

n  =1 


(4.6c) 


Z  A„{,)Q„(0)sinknx 


-  I  Al*  >  Wn  (0)  sin  kn  x  =  0. 


t  <  x  <  L 


(4.6d) 


where 


Pn  (0)  =  6r  — ~ —  coth  a„a)d  +  — — - , —  coth  an(2)  (h  —  d ) 

K  1-0'  t 


■02  1 
«i21  1  ~  02 


coth  aa(I)  {h  —  d  ) 


T„  {0)  =  0  2  — coth  a^’d  +  — Ifr  coth  a‘2)(/i  -  d  ) 


(4.7a) 


(4.7b) 
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Qn  W  = 


1  ~6r 

i -V 


coth  a„(2)  (h  —  d) 


(4.7c) 


(0)  = 


«r 


^  coth  +  -^y  coth  a„(2)  {h  —  d) 


1-02 


(4.7d) 


and  0  =  0/*o  .  This  set  of  equations  could  be  solved  in  their  current  form,  but  there  would 
be  no  savings  in  computation  time,  since  no  effort  has  yet  been  made  to  reduce  the  number 
of  terms  needed  for  the  summation.  Let  us  now  consider  the  method  for  doing  so. 

If  we  differentiate  Equation  4.6a  with  respect  to  x  and  substitute  the  result  into 
Equation  4.6c.  we  obtain 


£  Anif)kn  sin  knx  -  0.  0  <  x  <t 

n  =1 


£  A„ih  *  k„  sin  knx  =  0,  0  <  x  <  c 

n  =1 

Next  if  we  differentiate  (4.6d)  and  rearrange,  we  find 

£  A„(t )  kn  cos  kn  x  =  /  (x  ).  t  <  x  <  L 

n  =1 


£  A„(h )  ic„  cos  k„x  -  g(x).  t  <  X  <  L 

n 

where 


f  (x)  =  £  (am  Aj,‘ )  +  bm  Aj,h  >)  cos  km  x  . 

m  =1 


g(x)=  £  (cmAj,,)  +  dn,Aj,l,))  cos  kmx 

m  =1 

and 

Pn(0)W(0)-T(0)Q„(0) 
P(0)W(0)-T(0)Q(0) 


i 


(4.8a) 


(4.8b) 


(4.8c) 


(4.8d) 


(4.9a) 


(4.9b) 


(4.10a) 


i 


» 


W(p)=t - L.+1  (4. lid) 

1  ~*02 

It  turns  out  that  am.bm,cm.  and  dm  can  be  written  in  terms  of  expressions  such  as  P ~Pm 
and  T-Tm  .  and  that  these  expressions  decay  rapdily  for  large  m.  Thus,  for  example 

P  -  Pm  —  2er  +  c-<m-‘i)«A-rf)/£  (4.12) 

for  large  m.  Since  these  terms  decay  very  rapidly,  the  am  .  bm ,  cm  .  and  dm  coefficients  in 
Equation  4.9  also  decay  rapidly.  Therefore,  only  a  small  number  of  terms  are  necessary  for 
the  summations  in  Equation  4.8.  This  is  the  reason  that  the  method  is  numerically  very 
efficient. 

There  is  a  standard  method  of  solving  an  equation  of  the  form  of  Equation  4.8. 
which  appears  in  a  number  of  texts  [46-49],  as  well  as  in  [5].  We  therefore  forego  the 
detailed  explanation  of  the  solution,  and  instead  quote  the  result  from  [5].  The  solution  to 
Equation  4.8  is  the  following  matrix  equation. 


£(*,8,m  -  am  Dpm  —  Mm  Kp  )A„ 

m  *1 


-  L  (ft-  Dpn  +  Nn  Kp  )A,'A  >  =  0,  P  =  1.2. 


£  (— cm  Dqm  -  Xm  Kq  )Aj,'  >  (4. 1 3) 


+  Z  (  ft,  Dqn  -  Yn  Kq  )A„<*  >  =  0  ?  =  1.2,... 


where  8nm  is  the  Kronecker  delta.  As  before,  in  the  spectral  Galerkin  technique,  the 
solution  for  0  is  found  when  the  determinant  of  the  coefficient  matrix  vanishes. 

Let  us  identify  the  missing  terms  in  Equation  (4.13).  It  turns  out  that 
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^  ^mq  Iq 
q  *0 

z 


m— l 
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where  is  identified  by 


co  skmx  "-i 

- - -  =  L  pmq  C0S 

COS  kxx  9  =0 


(4.14a) 


(4.14b) 


(4.15) 


W-— g  ■ 

o  V 1  —  c«i  —  a2  cos0 


(4.16) 


y _ cos  9  d0 _ 

0  Vl  “  <*i  —  <*2  cos  ® 
r  sin  sin  0  d0 


q  =0 


q  -  1,2 . m  —1 


We  note,  furthermore,  from  Equation  (4.13)  that 
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V 1  —  oil  —  a2  cos  0*  —  Vl  _  <*i  ~ 


«2 


sin  ^0’  sin  0' d©’ 
Vl  ~  0[  -  a2  cos  0 


(4.18a) 

(4.18b) 

(4.19a) 

(4.19b) 

(4.19c) 

(4.20) 

(4.21) 

(4.22a) 

(4.22b) 


Finally,  the  last  two  terms  in  Equation  4.13  that  need  to  be  defined  are  K„  and  Dnm  .  These 
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2  Lr  cos  £  ,jc  sin  *„  x 

*■ =  x  •!  — a?? — dx 


Dnm  -  -i /cos*,x  L  f™,  sin  qQ  ~  Pmll 

**  I  q  =1  51 


sin  0 


(4.23) 


'  sin  kn  x  dx  (4.24) 


where  the  variable  0  is  defined  by  the  coordinate  transformation 


if  X  |  A 

cos  =  c«i  +  a2  cos  0 


(4.25) 


a,  =  Vz  cos  — -  —  1 


(4.26a) 


at2  =  ‘/z  cos  —  +  1 

Lt 


(4.26b) 


This  completes  the  description  of  the  final  solution,  which  was  stated  in  Equation  4.13. 

The  integrals  shown  in  Equations  4.16  and  4.17  are  well  behaved,  and  can  be  easily 
calculated  with  a  numerical  quadrature  subroutine.  The  integrals  in  Equations  4.20  and 

4.22.  however,  have  an  integrable  logarithmic  singularity  at  0  =0.  We  therefore  ^ust  use 
a  small  angle  approximation  for  cos  0  inside  the  logarithm  in  the  region  0  <  0"  <  9p . 
where  Qp  «  1.  and  use  an  integration  routine  in  the  region  0p  <  0"  <  tt. 

We  now  have  only  to  calculate  P ^  ,  Kn  ,  and  Dnm  .  as  shown  in  Equations  4.15. 

4.23.  and  4.24.  In  [5]  these  terms  were  calculated  only  to  the  second  or  third  order.  We 
now  would  like  to  extend  the  work  of  [5]  by  calculating  these  terms  out  to  the  fifth  order. 
These  are  calculated  in  the  next  section. 


4.3  Calculation  of  Pmq,  Kn.  and  DBm 

In  reference  [5],  Mittra  and  Itoh  calculated  P ^  .  K„  ,  and  D„„  to  the  second  or  third 


order.  We  would  like  to  extend  this  work  by  calculating  these  terms  to  the  fifth  order.  For 
convenience,  we  repeat  the  equations  from  the  last  section  in  which  these  quantities  were 
defined.  Hence. 


cos  km  x 
cos  k  \X 


=  £  Pm,  COS  96 
q  =0 


T  *  A 

2  r  cos  k  i*  sin  kn  x 
L  -  sin  9 


(4.15) 


(4.23) 


=  1/ 
L  J, 


COS  k  yX 


m  -l 


£  sin  qe  -  pm0  ^ 


sin  k.  x  dx 


(4.24) 


The  calculations  of  P and  Kn  are  straightforward,  although  somewhat  tedious.  The 
results  are 


*10-1 

P 20  ~  2a 2  1 

P2l  =  2o2 

P3  0  —  4af  —  2a!  —  1  +  2a22 
P 3i  =  -  2a2 

Pi2  =  2*i 

P40  —  8a]3  —  4a2  —  4a2  +  1  +  I2aja22  —  2a22 
P4l  ~  24a  2a2  —  8a2a2  —  4a2  +  6a23 
P4  2  —  12aia2  —  2a22 
?43  =  2a2 

jP50  =  16a!4  —  8a3  —  12a2  +  4ai  +  1  +  48a2af  —  I2axa$  —  6a|  +  6a24 
P  j]  —  64a3a2  —  24a  2a2  —  24aja2  +  4a2  +  48aja2  —  6a2 
Pi2  -  48a  2a22  —  12ata^  —  6a j  +  8a24 
P  53  =  16ata23  -  2a23 


K  i  -  a2 

K  2  —  <>2(1  +  2ckj) 

=  a2(4af  +  2aj  -  1  +  2aj 
K  4  =  a2(— 1  —  4a  3  —  4a2  +  2a2  +  8a2  +  120^  2) 

A'j  =  a2(l  —  4a  4  —  12a2  —  6a2  +  8a2  +  12a1a2  +  16a4  +  48a2a2  +  6a2)  (4.28) 

for  Kn  .  where  a3  and  a2  were  defined  in  Equation  4.26. 

The  calculation  of  Dnm  requires  a  rearrangement  of  terms,  because  the  algebra 
becomes  somewhat  involved  above  n-2.  Proceeding  from  Equation  4.24.  it  can  be  shown 


2a2  n— 1.2,... 

Dnm  ~  ——  L  L  CnpPmqlqp  m  -  12... .M 
/>30  ?  =0 


(4.29) 


where  Cnp  is  defined  by 


sin  k„  x  cos  k  ,x  " -i 

- =  £  CnP  cos'’  0  n  =  1.2,... 

sin^l  p  ~° 

Lp 


(4.30) 


and  Iqp  is  given  by 


j  —  sin  q9  cosP  *1  9  d9  q  —0 

o 

tr 

J  sin  q9  sin  9  cos^  9  d9  9  —1.2....  />— 0.1.2.... 


(4.31) 


The  elements  Cnp  and  lqp  are  considerably  easier  to  calculate  than  Dnm  .  We  obtain  for  C„ 
C  io  =  fa 


C,o  =  «,  +  ;/2 


C21  =  a2 


C3 o  =  <*i  -  &  +  2a22 

C  31  =  &2  +  40j(>2 

C  32  =  2a22 

C40  =  — %  —  2aj  +  2a2  +  4a)3 
C  41  —  — 2a2  +  4a)02  +  12aj 2a2 
C  42  ~  2 o<22  +  12a  )a| 

C  43  =  4a3 

C  so  =  1/2  —  2ai  —  6a)2  +  4a3  +  &a* 

CS1  -  —2  a2  —  12a)a2  +  12a)2a2  +  32a?a2 
C  j2  =  —6a2  +  12aja2  +  48a  2a22 
C  53  =  4a  3  +  32a)a3 


(4.32) 


and  the  values  of  Ifp  are  listed  in  Table  4.1.  By  using  this  method  of  generating  D„„  .  we 
reduce  slightly  the  efficiency  of  the  program,  by  introducing  a  double  summation  over  p 
and  q.  We  save,  however,  a  great  deal  of  algebra  that  would  otherwise  have  to  be  carried 
out  by  hand,  so  it  seems  like  a  small  price  to  pay. 

Mow  that  the  singular  integral  equation  technique  has  been  formulated  and 
extended  to  the  fifth  order;  we  consider  some  results  in  the  next  section. 


Table  4.1.  Values  of  /fP  as  defined  in  Equation  4.29. 
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4.4  Results 

The  propagation  constant  of  the  dominant  mode  of  a  shielded  microstrip,  as  shown 
in  Figure  4.1.  was  calculated  with  the  singular  integral  equation  technique.  A  comparison 
is  made  to  results  generated  from  the  spectral  Galerkin  technique,  and  to  the  results  from 
Mittra  and  Itoh  [5].  These  data  are  shown  in  Table  4.2  for  three  different  frequencies.  We 
can  make  several  observations  about  these  data. 

It  is  clear  that  there  is  very  good  agreement  between  the  calculations  using  the 
singular  integral  equation  technique  and  the  spectral  Galerkin  technique.  We  note, 
furthermore,  that  the  results  are  in  agreement  to  within  about  four  significant  figures.  This 
extends  somewhat  the  accuracy  achieved  by  Mittra  and  Itoh.  whose  results  are  shown  in 
Table  4.2.  This  increase  in  accuracy  is  due  to  the  larger  matrix  size  that  was  used  in  our 
calculations. 

The  excellent  correlation  between  the  results  generated  by  the  singular  integral 
equation  and  spectral  Galerkin  techniques  is  noteworthy  because  the  two  methods  are  of 


Table  4.2.  Calculation  of  the  dominant  mode  propagation  constant  or  shielded  microstrip 
as  shown  in  Figure  4.1.  These  calculations  were  made  with  the  singular  in¬ 
tegral  equation  (SIE)  technique  and  the  spectral  Galerkin  (SG)  technique.  The 
notation  2x2  indicates  the  matrix  size  in  the  SIE  technique,  while  the  notation 
2/50  indicates  that  two  basis  functions  and  50  spectral  terms  were  used  in  the 
SG  technique.  For  these  calculations,  h  -  2.0  mm.  L  -  1.75  mm.  t  «  d  -  0.5 
mm.  and  er  -  9.0. 


3  (rad/m)  at 

10  GHz  20  GHz  30GHz 

SIE. 

2x2 

525.574 

1094.59 

1685.50 

this 

4X4 

530.007 

1108.09 

1713.74 

paper 

6x6 

530.042 

1108.25 

1714.26 

10X10 

530.065 

1108.38 

1714.61 

SIE.  Mittra 

2x2 

530. 

1100. 

1710. 

and  Itoh  [5] 

4X4 

531. 

1115. 

1740. 

SG,  this 

2/50 

530.272 

1108.94 

1715.56 

paper 

5/100 

530.166 

1108.66 

1715.10 

I 
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very  different  natures.  It  was  not  clear,  based  on  previous  work,  that  the  spectral  Galerkin 
technique  could  be  used  to  generate  the  dominant  mode  propagation  constant  to  such  a  high 
degree  of  accuracy. 

An  attempt  was  made  to  calculate  higher-order  evanescent  modes  with  the  singular 
integral  equation  technique.  This,  in  theory,  is  easy  to  do  by  changing  the  propagation 
constant  from  purely  real  to  purely  imaginary.  In  practice,  it  turns  out  that  the 
determinant  function  that  must  go  through  zero  is  ill-behaved,  and  does  not  have  any  clear 
zeros.  Thus,  the  method  doesn't  seem  to  work  for  evanescent  modes.  One  explanatin  for 
this  result  could  be  that  we  have  only  calculated  enough  matrix  elements  to  fill  a  10x10 
matrix.  We  could,  of  course,  have  calculated  more  elements,  but  the  amount  of  algebra 
required  to  do  so  tends  to  become  very  large,  very  fast.  This  places  a  practical  limit  on  the 
matrix  size. 

4.5  Conclusion 

Although  the  original  goal  of  the  singular  integral  equation  technique,  that  of 
finding  higher-order  evanescent  modes,  was  not  met.  we  did  achieve  excellent  success  in 
verifying  the  propagation  constant  of  the  propagating  modes.  In  doing  so.  we  have 
demonstrated  agreement  between  two  very  different  methods,  the  singular  integral  equation 
method  and  the  spectral  Galerkin  technique,  to  four  significant  figures  in  the  final  answer. 
This  degree  of  accuracy  is  very  high,  given  the  complexity  of  the  calculation  and  the  very 
different  natures  of  the  two  methods  used  for  the  comparison. 
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CHAPTER  5 

COPULED-MODE  ANALYSIS  OF  MULTICONDUCTOR  MICROSTRIP  LINES 

5.1  Introduction 

Recent  advances  in  microelectronic  packaging  have  generated  a  number  of  difficulties 
associated  with  interconnections  between  VLSI  logic  devices.  These  interconnections  are 
typically  made  with  large  numbers  of  microstrip  lines  that  run  parallel  to  each  other.  It  is 
possible,  under  certain  conditions,  that  a  signal  propagating  on  one  line  can  couple  to  other 
lines.  This  cross-talk  can  trigger  false  signals  on  nearby  lines.  Therefore,  it  is  important  to 
fully  analyze  multiple  microstrip  lines,  in  order  to  determine  the  performance  of  the 
overall  system.  These  cross-talk  problems  become  more  serious  as  switching  speeds  and 
packing  densities  are  increased. 

Coupled  microstrip  lines  have  been  analyzed  with  several  different  methods 
[6.7.50.51].  More  recently,  multiconductor  transmission  lines  with  more  than  two  lines 
have  received  some  attention  [52.53],  These  analyses,  however,  are  all  based  on  a  quasi¬ 
static  approximation  which  is  valid  only  for  digital  devices  with  switching  speeds  in  the 
order  of  one  nanosecond.  When  the  rise  time  of  the  switching  pulse  is  reduced  to  hundreds 
of  picoseconds,  the  spectral  content  of  the  pulse  contains  higher-frequency  components. 
Since  quasi-static  approximations  break  down  at  high  frequencies,  a  full-wave  analysis  of 
multiconductor  microstrip  transmission  lines  becomes  necessary. 

The  configuration  under  study  is  shown  in  Figure  5  1  It  is  composed  of  five 
microstrip  lines,  each  of  the  same  width  and  strip  separation.  In  this  chapter  we  would  like 
to  demonstrate,  using  a  frequency-dependent  analysis,  the  mechanism  whereby  current  that 
is  excited  on  one  line  is  transferred  to  other  lines.  We  begin  with  an  analysis  of  the  modes 
of  propagation  on  the  line. 
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5.2  Calculation  of  Modes 


In  a  multiple  microstrip  configuration  with  N  strips,  there  are  N  modes  of 
propagation  down  these  lines.  In  order  to  analyze  the  crosstalk  and  propagation  delay  of 
these  lines,  the  propagation  constant  for  each  of  the  N  propagating  modes  is  needed.  The 
method  we  use  for  this  calculation  is  the  spectral  Galerkin  technique,  described  in  Chapter 
2.  By  making  a  number  of  alterations,  the  same  theory  can  handle  multiple  lines  as  well  as 
a  single  one. 

We  begin  the  analysis  with  the  matrix  equation  generated  by  the  spectral  Galerkin 
technique  for  a  single  shielded  line.  This  equation  was  given  previously  in  Equation  2.16. 
for  the  configuration  shown  in  Figure  2.1.  In  order  to  remove  the  top  wall,  we  make  a 
modification  in  the  spectral  Green’s  function.  This  is  done  by  replacing  the  short  circuits  in 
the  equivalent  transmission  lines,  at  y  -  h  -  d.  shown  in  Figure  2.3.  with  matched  loads. 
The  side  walls,  on  the  other  hand,  are  removed  by  moving  them  far  enough  away  from  the 
strips  that  they  have  no  effect.  Care  must  be  taken,  however,  not  to  move  them  too  far 
away,  since  the  number  of  spectral  terms  required  for  this  calculation  can  become  large 
with  wide  side  walls. 

All  that  remains  now  is  to  generate  a  reasonable  set  of  basis  functions  for  the 
current  on  the  lines.  We  choose  to  build  our  basis  functions  up  from  a  set  of  functions  that 
would  be  suitable  for  finding  the  odd  and  even  modes  in  a  microstrip  containing  only  a 
single  strip.  These  building  block  functions  are 

e  (,  \  _  cos  ((i  -  1)tt(x  Is  +  1))  _  ^  _  sin  ((t  -  Vzfrix  /s  +  1)) 

'•**  /  "  1  -  ifet  \X  )  —  - 

v  1  -  (jc/j  )2  V 1  -(x  Is)2 

(5.1) 

e  ,  _  cos  ((t  —  Vz)ir{x  Is  +1))  .  .  _  sin  (i  n(x  /s  +  1 )) 

So/  ' *  ’  ,  = -  ~oi  i  *  '  —  -  - — 

V 1  —  x  Is  )2  vl-U/t)2 

We  note  that  and  r)0,  are  the  same  basis  functions  used  for  the  even  modes  in  a  single 
line  in  Chapter  2.  These  were  given  in  Equation  2.22  and  shown  in  Figures  2.6  and  2.7. 

We  now  take  advantage  of  the  symmetry  inherent  in  the  structure  to  find  basis 


functions  that  are  either  even  in  Jz  and  odd  in  Jx  .  or  odd  in  Jz  and  even  in  Jx.  Although 
the  method  that  follows  is  general  for  N  lines,  we  show  as  an  example  the  basis  functions 
used  for  five  lines  as 


V:  5 

/*(*)  =  £  £  c,j  ptJ(x ) 
;=li=l 

5 

A  (* )  =  £  £  d,j  Qij  (x  ) 

y=l/=l 

where,  for  even  Jz  and  odd  Jx 


P  !,•(*)*  t,  (*  > 

/>2l  (x  )  =  £tl  (x  +  q  )  +  i„  (x  -q) 
Pi,  (*  )  =  L,  (x  +  q)~L,(x  -q) 
Pa,  (*  +  2q)  +  £„  (x  -2 q) 

P$,  (x  )  =  |01  (x  +2  q)-  £ol  (x  -  2?  ) 

and  for  odd  Jz  and  even  Jx  . 

Pi,  (x  )  =  £01  (x  ) 

Pi,  (x  )  =  t,  (■*•  +  <7  )  “  £«  (*  -  <? ) 
Pi,  (x  )  =  £0,  (x  +q)  +  f0,(x  -  q) 
Pa,  ( x  )  =  (x  +  2?  )  -  £„  (x  -  2q  ) 
P5l  (x  )  =  (x  +  2q  )  +  £01  (x  -  2q  ) 

where 

<7  =  2  (i  +  r  !) 


Q  li  (x  )  T)0^  (x  ) 

Qi,  (x  )  =  Vo,  (x  +q)  +  Vo,  (x  ~  q  ) 
Qi,  (x  )  =  Vt,  (x  +q)-r (x  -  ?  ) 
Q4,  (x  )  =  t)oi  (x  +  2?  )  +  Vo,  (*  ~  2q  ) 
<2  5,  (x  )  =  Vt,  (*  +  2q  )  -  Vr,  (x  -2q) 


Q  i,(*)  =  Vf,  (x  ) 

Q  2,  (x  )  =  1%,  (x  +  )  -  T70,  (x  -  9  ) 

£?3,  (x)-vf,(x  +  q)  +  TJ„  (x  -  ?  ) 

(2  4,  (x  )  =  1%,  (x  +  2if  )  -  Vo,  (x  -2 q) 

Q  s.  (x  )  =  T)„  (x  +  2?  )  +  7},,  (x  -  2?  ) 


(5.2) 


(5.3a) 


(5.3b) 


(5.4) 


With  these  basis  functions,  we  now  search  for  the  three  even-mode  and  two  odd-mode 
propagation  constants  that  satisfy  Equation  2.16.  Results  showing  dispersion  curves  and 
comparisons  to  the  quasi-static  theory  are  presented  in  Section  5.4. 
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Let  us  now  turn  to  the  coupled  mode  theory  used  to  explain  the  transfer  of  current 
between  lines. 

S3  Coupled-mode  Theory 

When  a  single  line  of  an  N'-line  system  is  excited,  one  may  consider  this  to  be  an 
excitation  of  a  linear  combination  of  the  five  modes,  whose  weighting  coefficients  must  be 
determined.  We  may  consider  the  integral  of  Jz  over  an  individual  strip  to  be  the  quantity 
that  is  coupled  between  the  lines.  Therefore,  we  carry  out  a  procedure  similar  to  that 
described  by  Bhartia  and  Bahl  [54]  for  image  guide  couplers.  For  each  mode,  we  find  the 
longitudinal  current  on  strip  i  due  to  mode  j  as 

fj,  =  f  Jt(x)dx  i.j  =  1 . 5  (55) 

strip  t 

We  consider  an  excitation  of  the  first  line  alone  to  be  one  where  the  first  line  has  a 
magnitude  of  one  and  the  others  have  a  magnitude  of  zero  at  z  -  0.  Therefore,  we  set  up  a 
fifth-order  linear  equation  to  find  the  combination  of  modes  that  satisfies  the  condition 

L  aiIP  = 

i  =  i 

where  the  a,  's  are  the  unknown  mode  coefficients. 

Once  the  mode  coefficients  are  determined,  we  may  assume  they  propagate  down  the 
five-line  system,  with  their  respective  propagation  constants  and  mode  coefficients  that  were 
calculated  in  the  previous  section.  Thus. 

/;(*>«  I  a,  /„  e (5.7) 

where  is  the  propagation  constant  of  the  j'h  mode,  and  I‘(z  )  is  the  current  on  the  i'h 
strip  due  to  all  j  modes,  at  a  distance  z  from  the  original  excitation.  This  can  be  calculated 
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for  each  strip  as  a  function  of  2.  and  results  are  shown  in  the  next  section. 


S.4  Results 

Let  us  now  examine  the  results  obtained  with  the  method  described  above.  The  first 
step  of  the  procedure  involved  calculating  the  propagation  constants  of  the  various 
propagating  modes  associated  with  the  configuration  in  Figure  5.1.  We  have  calculated  the 
modes  at  1  GH2  for  a  three-line  configuration  using  the  spectral  Galerkin  technique,  and 
then  we  compare  them  to  those  of  Chan's  quasi-static  approach  [53],  This  comparison  is 
shown  in  Table  5.1.  At  this  comparatively  low  frequency,  the  two  methods  are  in  excellent 
agreement.  At  higher  frequencies,  however,  the  quasi-static  method  is  expected  to  break 
down. 

This  point  is  demonstrated  in  Figure  5.2.  where  the  five  propagation  constants  of  a 
five-strip  system  are  plotted  as  a  function  of  frequency.  At  low  frequencies,  the  curves  are 
level,  as  would  be  expected  by  quasi-static  theory.  At  higher  frequencies,  however,  the 
relative  propagation  constants  squared  are  no  longer  constant,  indicating  that  a  frequency- 
dependent  theory  is  now  necessary. 

Next,  we  checked  the  convergence  of  the  calculated  propagation  constants  with 
respect  to  the  number  of  modes.  These  results,  shown  in  Table  5.2.  demonstrate  a  very 
rapid  convergence  with  respect  to  the  number  of  modes  and  suggest  that  probably  one  basis 
function  of  each  type  (A/f  =  iV,  =1)  will  be  sufficient  for  most  calculations. 


Table  5.1.  Comparison  of  propagation  constants  calculated  by  two  different  methods. 

Shown  here  are  the  two  even  and  one  odd  modes  of  a  three-microstrip  system 
at  1  GHz,  where  t  -  2s  -  1  mm.  and  2s  1  -  0.2  mm. 


Spectral  Galerkin 
0,  i(rad/m)  50.53 

0,2  58.12 

0„  52.80 


Juasi-static 
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Dispersion  curves  for  the  first  five  inodes  of  propagation  of  a  microstrip  with 
five  lines.  The  dimensions  are  the  same  as  those  in  Figure  5.1. 


Table  5.2.  Variation  of  the  five  propagation  constants  in  a  five  line  system  with  the 
number  of  basis  functions.  For  this  configuration,  t  *  2s  -  1  mm.  2rj  -  0.2 
mm,  and  f-1  GHz. 


EKEfl 

ft. 

flUIKi 

■■ 

i 

49.86 

50.83 

60.50 

49.97 

53.58 

2 

49.45 

50.77 

60.40 

49.75 

53.57 

3 

49.42 

50.77 

60.40 

49.73 

53.57 

In  order  to  get  a  physical  feel  for  the  shape  of  the  five  modes,  we  have  plotted  Jz  as 
a  function  of  x  over  a  cross  section  of  the  strip  for  each  of  the  five  modes.  These  are  shown 
in  Figures  5. 3-5. 7  for  Nz  —  N,  =1  and  for  Nz  —  Nx  —  3.  We  note  that  for  most  of  the 
modes,  we  may  obtain  a  very  reasonable  representation  of  the  currents  with  just  one  basis 
f unction  of  each  type. 

Finally,  we  demonstrate  the  coupling  of  current  from  one  strip  to  the  next.  In 
Figure  5.8.  we  present  results  for  the  configuration  of  Figure  5.1  at  1  GHz.  and  of  Figure 
5.9  at  10  GHz.  We  start  with  a  unit  excitation  on  the  first  strip  and  note  that  the  current 
on  the  first  strip  very  rapidly  decreases,  while  the  current  on  the  other  strips  rapidly 
increases.  Thus,  when  microstrips  are  closely  spaced,  the  coupling  between  strips  is 
predicted  to  be  quite  severe. 


5-5  Conclusion 

In  this  chapter,  a  frequency-dependent  method  of  calculating  the  coupling  between 
a  large  number  of  parallel  microstrips  has  been  demonstrated.  Using  this  method,  it  has 
been  shown  that  coupling  between  closely  spaced  lines  can  be  quite  severe.  This  has 
significance  in  the  area  of  VLSI  interconnectivity,  where  parallel  microstrip  lines  are  packed 
as  densely  as  possible. 


J,(x) 


figure  5.8. 


z  in  free-space  wavelengths 

Variation  of  /  with  z  at  1  GHz  for  each  of  the  five  strips. 
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Variation  o!  /  with  z  at  10  GHz  for  each  ol  the  five  strips 


CHAPTER  6 


SUMMARIES  OF  OTHER  ACTIVITIES  AND  IMPORTANT  RESULTS 

Other  activities  in  this  effort  that  have  been  reported  previously  in  various  technical 
reports  have  included  the  study  of  coplanar  transmission  lines  [55.  56]  and  fin  lines  [57.  58] 
and  the  investigation  of  dielectric  antennas  for  millimeter-wave  imaging  applications  [59]. 
In  addition,  a  study  of  planar  waveguides  and  components  for  millimeter-wave  integrated 
circuits  has  been  carried  out  [60].  In  Ref.  60,  the  problem  of  discontinuities  in  planar 
guides  has  been  treated  using  the  generalized  variational  method  and  the  conjugate  gradient 
technique. 

Summarizing  the  important  results,  we  have  shown  that  the  spectral  Galerkin 
method  is  useful  for  calculating  the  evanescent  as  well  as  propagating  modes  in  a  shielded 
microstrip.  We  have  demonstrated  that  the  accuracy  of  these  modes  is  quite  reasonable:  the 
field  distributions  corresponding  to  these  modes  have  been  plotted  to  provide  a  physical 
insight  of  their  structures. 

Next,  these  modes  have  been  used  in  a  mode-matching  analysis  to  calculate  the 
scattering  from  discontinuities  in  various  microstrip  configurations.  Again,  reasonable 
results  have  been  obtained  for  single  and  cascaded  discontinuities:  however,  the  problem  of 
accurately  calculating  the  junction  capacitance  has  presented  some  difficulties  that  have  yet 
to  be  overcome. 

We  conjecture  that  this  is  due  primarily  to  our  inability  to  calculate  a  sufficient 
number  of  modes  with  enough  accuracy,  as  evidenced  by  our  check  on  the  orthogonality  of 
these  modes.  In  addition,  the  matrix  generated  by  these  modes  was  also  found  to  be  ill- 
conditioned. 

Next,  we  investigated  an  alternate  method,  viz.,  the  Singular  Integral  Equation 
technique,  for  calculating  the  modes  in  a  shielded  microstrip  line.  The  calculated  values  of 
the  propagation  constant  of  the  dominant  mode  were  in  agreement  to  four  significant 


figures,  with  results  generated  with  the  spectral  Galerkin  technique.  This  degree  of 
accuracy  is  exceptional,  because  the  two  methods  that  are  being  compared  are  quite 
different.  An  attempt  was  made  to  calculate  evanescent  modes  using  the  singular  integral 
equation  technique,  but  these  modes  were  not  found.  One  possible  explanation  was  that  the 
size  of  the  matrix  we  could  use  was  too  small.  Using  a  larger  matrix  was  not  practical 
because  the  complexity  of  the  matrix  elements  increases  rapidly  with  matrix  size. 

The  generalized  variational  method  and  the  conjugate  gradient  method,  both  of 
which  are  iterative  techniques,  have  been  employed  to  derive  the  solution  to  the 
discontinuity  problems  in  planar  guides,  but  they  have  also  yielded  only  limited  success  in 
terms  of  the  realized  accuracy  of  the  results.  Thus,  our  conclusion  is  that  further  effort  in 
this  direction  is  critically  needed. 

Finally,  the  dielectric  antenna  we  have  developed  has  been  demonstrated  to  be 
useful  as  an  array  element  of  an  imaging  system  at  80  GHz  and  220  GHz.  Again,  a 
monolithic  imaging  system  would  be  better  suited  for  this  purpose  and  it  would  be 
worthwhile,  in  the  future,  to  investigate  such  a  monolithic  system. 
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